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ATTENTION 

An error was made in not numbering pages in this booklet. A few moments on the 
part of the user will take care of this oversight with use of this Table of Contents 
as a guide. 



The purpose of this pamphlet is to present, in a consistent form, the Statisti- 
cal Library Routines available to Digital Computer Laboratory users on Illiac. 
The condensed routine writeups -will give, when available, information ons 

lo General purpose of routine . 

2. Operator's Instructions (This portion may be copied on Operator's 

Instruction sheet submitted with production problem). See (Procedure 

for submitting a production problem). 
3» The preparation of tapes to be used. 
ko Estimated time for running problems 
5o FF stops (usually indicative of an error) if any. (This should also 

be copied on Operator's Instruction sheet). 
6. Routines, if any, that will utilize the output with little or no 

change on output tape. 
7- Gomments that, it is hoped, will facilitate the use of these routines. 

In order to gain consistency throughout, the following nomenclature will be used. 

R-Matrix s An R-matrix is a correlation matrix obtained generally from programs 
labeled K-2, K-8, or K-9. It can be obtained in two forms; the lower triangular 
or in rectangular form. IMPORTANT : If the intention is to use the R-matrix for 
further Illiac use e.g. v-f actor: analysis,^ it: -must be in triangular form. 

ILLIAC Row ; Most Statistical Routines use a form of Illiac readin that places 
an imaginary decimal point immediately following the sign. This should be held 
in mind when scaling. For example, K-8 output is scaled 10" 1 . The output tape 
may be — - +1000 -0293 +0026. This would be interpreted as +lj -.293; +.026. 

It is possible to place a group of numbers on Illiac tape such that your 

teletype printup will be in columnar +029 or row, +029 -032 +046, form. (Or 

-032 

+046 
conceivably, scattered over the paper +029 n ^ Q ). To Illiac, these numbers 

" U ^ +0k6 
are considered as a row of numbers when it is followed by an N, J, F, or L, 

^th-Hole Gharacter ; At times a fifth-hole character (i.e., carriage return- 
line feed 1 space j one-hole delay) is an integral part of a parameter or data tape. 
When such a character is required it will be so indicated, and must be included. 
If no indication is made, these fifth-hole characters can be included or may be 
excluded without resultant harm. 

+ or Kg - or S g Unless specifically indicated, a + and - may be punched 
instead of the letters K and S (e.g., a parameter requiring a no^ S may be punched 
as no. -)o 

n and v ; Much statistical work involves the use of a set of observations. 
The set may be people, observations or any group which contains for each unit in 
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that group a similar number of measurements or observations for each sample « 
These are also called variables or tests. Unless otherwise noted in this pamphlet 
the measurements j regardless of what they are, will be referred to as v or variables 
and the group which provides the measurements will be referred to as n or sample 
size. In general, this will take the form of a rectangular matrix with n rows, each 
row consisting of v variables. At times the term f factors will be used. This 
consists of f columns with v variables in each column. 

d~ and d,. ; In specifying d . (decimal digit input or the number of digits in 

each element input), leading zeros as well as trailing zeros are to be counted. 
A number +00 300 would have an i_ value of 5« This holds true for dL (decimal digit 
output or the number of digits in each element output). Most- of these routines are 
scaled 10 . This means that a leading zero will be punched out. Therefore, if 
you wish 5 significant figures in your output, it is a safe thing to ask for a cL 
value of k. When exceptions to this arise, the program writeup will indicate them. 

Output Read Back in ILLIAC ; Many times you may wish to use Illiac output for 
input in a different problem at a different time. Be sure to check the output tape 
to be sure extraneous titles, headings, etc., are not present. Generally not doing 
this will cause Illiac to hangup. Either remove the offending part, or mark the 
tape where the operator should begin reading the tape. When practicable in the 
following descriptions, such output (with headings, etc.) will be indicated. 

Factor ; A column or row of _v_ (variables) signed fractions. 

R-Matrix ; Any correlation matrix of v (variables) by v signed fractions. 

Matrix ; An array of numbers placed in an n row, in column form. Usually 
denoted by A. ., where i = number of rows and j = number of columns. 

Output ; Results received from Illiac. 

Input ; Information on tape fed into Illiac. 

Triangular Matrix Form; The lower left hand side or upper right hand side 
of a symmetric matrix. 

Square Matrix Form ; An array of elements in matrix form where the number of 
rows and columns are equal. 

Routine ; Any tape that has been written to solve a problem on Illiac . 

Master Tape s Any tape that has been written to solve a problem on Illiac. 

Program Tape ; Any tape that has been written to solve a problem on Illiac. 

Variable ; A single measurement that can vary in value for each time measured. 

Sample Size; The number of times an identical set of variables are measured. 
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Order of Matrix ; . A square matrix of n rows and n columns has an order n. 

Columnar Form Printout : Results printed as a column on the teletype printer. 

Row Form Printout : Results printed as a row on the teletype printer. 

Matrix Form ; Results printed out as n rows and m columns in a rectangular form* 

Parameters : Used to limit Illiac programs to your particular problem. 

Signed Numbers ; Integers or fractions preceded by a + or -. 

d n : The number of digits (does not include sign) in each element of output - 
include all zeros . 

■d.: The number of digits (does not include sign) in each element of input - 
include all zeros. 

A. .; Used to denote a matrix with i rows and j columns. 



Black Stop ; Any 2.K, 20/$* 50 order. Illiac stops on these. * oh.ca,: black switch 
up causes Illiac to resume. 

Black Down ; Called "black disable", "black ignore", "stop disable" ^.causes 
Illiac to ignore or bypass all black stops. 

Terminating Symbol ; Any symbol (usually N, J, F, L) which ends a row or 
column of numbers and/or a matrix. 

Element : A signed fraction or integer, usually a. . in a matrix. 

Observations : Sample size, a set of variables measured over a group n of people. 

Tests : Variables or measurements. 

. 5th~hole Delay : An Illiac punch which has one hole in the fifth position. 
Does not show up on teletype printer. 

5th-hole Character : Any Illiac character other than thru L that contains 
a fifth-hole punch. 

Sexadecimal Character : The Illiac characters 0, 1, 2, 3, ^, 5> 6, 7, 8, 9, 
K, S, N, J, F, L. 

Character : Any column of i:ey punched holes on Illiac tape. 

Diagonal Element : In a matrix of n rows, n columns the elements a.^, a^ 3 
..., a , i.e., element at intersection of row 1, column .1; row 2, column 2, etc. 

Cumulative : Summed or accumulated values over calculations. 

Machine Error : Error on Illiac caused by machine malfunction. 

Programmer Error ; Error on Illiac caused by user. 

Absolute Value : Any number, positive or minus, that is assigned a positive 
value . 

Symmetric Matrix -;- A matrix A. . whose a . element = a element. 

— — — — — — — — — — — — — lj 1 J J -L 
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Transformation or Transform Matrix: A matrix x.. which transforms a matrix 
— ■ ■ Ji 

A., into matrix B. n • 
lk jk 

Rectangular Matrix ; A matrix of n rows and m columns . 

Vector : One row or column of a matrix. 

Fixed Number : Any number where the decimal point is fixed by user<> 

Operator's Instructions : Instructions filled out on an Illiac "Program 
Instruction" sheet .. to enable operator to run problem. 

Sums of Sows and Columns : Zr and Zc. A value obtained by summing elements 
across rows or down columns • 

Matrix Transpose s Take a matrix A. . and interchange rows and columns to get A... 

Scatter Plot : A graphic plotting of two columns of numbers ^ one plotted as 
X the other as Y on an XY axis. 

OF Stop : A stop on Illiac that usually designates the end of a machine run. 

Identity Matrix : : A matrix with ones in the diagonal elements and zeros for 
all other elements . 

Scaler Matrix: A matrix with the same constant in each diagonal element and 



zeros in all other elements . 



A. .+B. . 



a ll a 12 a 13 
a 21 a 22 a 23 
a 31 a 32 a 33 



b ll b 12 b 13 
b 21 b 22 b 23 
V b 32 b 33 



a li± b ll a 12± b 12 a 13 ±b 13 
a 21± b 21 a 22i b 22 a 2^23 
a 31 ±b 31 a 32 ±b 32 ^33^33 



Frequency : The number of times a number or. event occurs. 

Non-singular Matrix : A matrix whose determinant value = : 0. 

Capacity of Program : The highest values of data size that can be run on 
any one program. Do not exceed. 

- Sums of Squares of Rows and Columns : Zr and Zc . Each element in a row 
or cxiLuaan is squared. The squared elements are then summed over rows or columns 



A. .xB.. 
ij jk 



a ll a 12 a l3 
a 21 a 22 a 23 
a 31 a 32 a 33 



b ll b 12 b 13 
b 21 b 22 b 23 
b 31 b 32 b 33 



( a ll b ll +a 12 b 21 +a 13 b 3i } (a ll b 12 +a 12 b 22 +a 13 b 32 } (a ll b 13 +a 12 b 23 +a 13 b 33 } 
(a 21 b ll +a 22 b 21 +a 23 b 31 ) (a 21 b 12 +a 22 b 22 +a 23 b 32 3 (a 21 b 13 +a 22 b 23 +a 23 b 33 ) 
(a 31 b H +a 3^1 +a 33 b 31 ) (a 31 b 12 +a 32 b 22* a 33 b 32 ) (a 3 l b 13 +a 32^3^33^ 3 ) 
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- FFOOXX Stops : These stops ¥-111 indicate two types of errors to the user. 
They are; (l) a machine error or; (2) an error on the part of the user, i.e., 
data tape error, parameter error, or the structure of the data causes an arithmetic 
error. Should such an error occur, the operator will note the error on the 
'^Operator's Instruction Sheet" and return it to you. The type of error is designated 
by the XX and/ or location of the FF stop*- It is requested that you list, when 
given, on your "Operator's Instruction Sheet" the FF stops given in the following 
writeups. Also include a description, when given, of the errors indicated "by these 
stops o 

Time Estimates : All production problems submitted must have time estimates. 
These are used by the computer operator's in scheduling production. Since an 
underestimate of time will result in your problem being returned before it reaches 
completion, it is to your advantage to see that the time estimate is great enough 
to cover the running of your problem. Time estimates are given for routines when 
possible, -Many times these are, at best, estimates* 

Procedure on Submitting Production ; Following is a step by step procedure 
for using Illiac facilities: 

1. Fill out and submit to the- Illiac operator a "Problem Specification" 
sheet (this will be returned in a few days, and you will be assigned 

a problem number for use on Illiac)* Obtain the number from the operator „ 

2. Place data on tape according to the format required for the routine 
you are using* 

3. Place on tape any parameters that are required. 

k« Prepare the master routine for your problem. There are three 
libraries which contain copies of master routines. One is the 
Statistical Library (denoted by KSL); the Illiac Active Library 
(denoted by an alphabetic letter and number such as K-8, M-12, 
etc) and the Illiac Auxiliary Library, The Illiac teletype operator, 
can help you locate these. The tape in the library must be duplicated 
by you. This duplicated routine is then submitted with your production 
problem , 

5. Fill out a "Operator * s Instructions " sheet for the Illiac *: cAllc blank 
spaces are important- — -fill' them:' in « It is on these sheets that 

is placed the sections in these writeups called "Operator's Instructions" 
and "FF Stops". 

6, The operator will at this point give you the necessary information 
on: how to turn in your problem; when you can possibly expect your 
results; and where you can pick up your results. 

Black Switch, White Switch : On Illiac are two switches used by the operator 

in running a program. Generally, unless told differently on your ^Operator's 

Instructions" , the operator will use the white switch to read in the master and will 

use the black switch from then on. There are times you need to give further :' 

- 5 - 



instructions to the operator . If the writeup herein includes instructions to run 
black down or to white switch up and down to read in a new tape, this instruction 
is important, and must he included on your instructions-, The operator's are not 
mind readers <> Therefore, include any unusual instructions if you wish efficient 
service. If no instruction is given, the operator will assume you wish to use 
the black switch once your master tape is in. 

Data Input ; Many times a program will not read in all of the data tape. 
This occurs when output from one program is read back into Illiac as input on 
another program and only part of the output is needed for input • It is important 
that you indicate this to the operator. -Such indication might be 

(3) Data tape (will not all read in) 2k--- 

NOTE; AT TIMES, ERRORS ARE FOUND IN THE DESCRIPTIONS OF ILLIAG PROGRAMS. 
SHOULD YOU, IN USING ANY OF THE FOLLOWING ROUTINES, FIND ERRORS, PLEASE INDICATE 
THEM TO THE SECRETARY IN THE TELETYPE ROOM SO THAT SHE MAY MAKE THE NECESSARY 
CORRECTIONS . 



KSL 1.11 
PRINCIPAL AXIS FACTORS 



PURPOSE 



This routine will extract the largest K factors from a 
correlation matrix (K=2, K-8, K=9 output in triangular 
form) by an iterative procedure suggested by Hottellingo 
The routine will fail in the exceptional circumstance 
that two or more of the largest K eigenvalues are of the 
same size<> The number of variables in the correlation 
matrix must be < 111„ The output is in the form of a 
Title (must be eliminated for further runs of output 
on Illiac) the factors in columnar form, signed, scaled 
10°° 9 an N at end of each column <> A variance table giving 
Factor number , latent root or eigenvalue, percent of 
variance extracted for each factor, and the cumulative 
percentage of variance extracted,, This table is followed 
by the residuals left in communalities and a triangular 
matrix of residual variance left in original correlation 
matrix „ 



OPERATOR'S INSTRUCTIONS 



lo Master tape 

2. Parameter tape 

3» Communality tape 

h* Data tape (black switch down) 

5« (Center black switch during residual 
output ) 




Black switch up 



PREPARATION OF TAPES 



3^122 

2*LLN0 (if no 
residuals desired, 
stop at this point) 

At this point a new parameter tape, communality tape and 

data may be read in a 

lo Master - Duplicate library copy KSL loll 

2 o Parameter 

where F = number of factors desired 



F space d» space d space 



d n = number of decimal places in output 

d = number of decimal places in residual output. 



(This 



may be set to if no residual matrix is desired) 



TIME ESTIMATE 



FF STOPS 
OUTPUT 



3» Communality tape - this consists of a column of signed 

elements scaled 10 and terminated by an N. The elements 
may be either +1 or trial values less than +1. For 
values less than .1, a zero will precede each value. 

ko Data tape— this usually will be K-8 output. It consists 
of a triangular matrix,, each element scaled 10 and 
signed. If the matrix is punched by hand it consists 
of the upper half of a symmetric matrix by columns or 
the lower half by rows . 

No time estimate given- -varies greatly. 

For example, in a 5 factor, lk variable problem the following 

times were noted. 

Fl - 25 seconds^ F2 ~ 25 seconds! F3 - 85 seconds j.f4 - 30 

seconds ! F5 — k^ seconds— —30 seconds to punch out k decimal 

residual matrix. 

Two 50 variable 2 factor problems were run. Their times were z 

1) First factor - 1 minute 6 seconds 5 second factor - 1 minute 
^0 seconds and 3 minutes to punch out residual matrix. 

2) First factor - 50 seconds! second factor - 3 minutes 
hk seconds .- 3 minutes to punch out residual matrix. 

None listed 

The output consists of a Title! ■■"the factors printed in 
columnar form, sealed 10 ~ , each column terminated by an N| 
a variance table! an N terminated column of communality 
residuals! and a triangular (lower half) matrix of residuals. 
Before using the factors for a run in Illiac, be sure the 
title has been removed. 
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PURPOSE 



KSL 1.12 
M-22 OUTPUT TO FACTOR LOADINGS 

This routine takes M-22 output and converts it to the 
desired number of factors where the number of factors 
is equal to or less than the number of factors in M-22 
output. Output from 1.12 is scaled 10" , by columns , 
and has an N at the end of each column and a J after the 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



ESTIMATED TIME 



EXAMPLE 



last N. 

1. Master KSL 1.12 

2 • Parameter 

3. Data (M-22 Output) 



Stops 
2405J 
24078 
24078 



At stop 24078 raise black switch to read an additional 

data tape using the same parameters . 

To read a new set of parameters at 24 3 move the white switch 

up and down. 

Master - Duplicate Library Master Copy KSL 1.12 

Parameter - 003+ 

00F OOf F 

OOF OOvF 

OOF 00d Q F 

24120N 
where 
f = number of factors to be converted (equal to or less than 

number of factors in M-22 output) 
v = number of variables in each factor 
d = number of decimal digits in the output 
Data - Output from M-22 
Master - 1 minute to read in 
Data - approximately 150 characters/second 
Output - approximately 50 characters/second 
Data tape = 10 factors % 50 variables; 7 decimal places 

Data read in = 10 g0*7 ~ 2 6 seconds 



Data output = 10 factors; 50 variables; 10 decimal output 

Data output time = J? = 100 seconds 

Total time = 26 seconds + 60 seconds + 100 seconds = 3 
minutes* 
FF STOPS None 
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KSL 1.20 - 251 
CENTROID FACTORS WITH FIXED COMMUNALITIES 

PURPOSE This routine takes an R-matrix of v variables , (output 

from K-2, K-8, K-9) a set of communalities provided "by 
the user (of v size) and reduce the R-matrix into a set 
of f linearly dependent factors. The communalities in- 
troduced replace the diagonals in the R (correlation) 
matrix.* There is no re-estimation of communalities as 
in KSL 1.21. The resultant factor matrix satisfies the 

equation R = FF' + R, ^,\ where R/ »\ is the matrix of 
^ v (v-f) (v-f) 

residuals . This program is limited to 111 variables and 
99 factors. Except for a title at the front, results from 
this program may be fed directly into KSL 1.90, 1.80 with- 
out any modifications of the tape. The output contains 
a title, the factors signed, scaled 10 by column - an N 
at end of each column. The factors are followed by a 
variance table containing the factor number! the latent 
root or eigenvalue) the percentage of variance extracted 
for each factor and the cumulative percentage of variance 
extracted. This, in turn, is followed by a column of 

residuals left in the communalities and the R/ „\ matrix 

(v-f) 

which can be fed in later for further factor extraction. 

If this is done, the operator tape at the end of the output 
is necessary.** 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL 1.20 2^190 

2. Parameter tape 2^l4L 
5 • Communality tape 3^1 $5 
k. Data tape (correlations) 2^190 

A new parameter tape.,_ communality tape, and R-matrix 
may now be read in without rereading the master. 
TAPE PREPARATION 1. Master tape - Duplicate library copy KSL 1.20 

2. Parameter tape - this consists of three unsigned decimal 
numbers separated by a fifth-hole character (space, 



carriage return, delay) they are f space d space 

d space. - 

r 

Where f is the number of factors desired, (must be < 100) 

d„ is the number of decimal places wished in each 

element of the factors and 

d is the number of decimal places wished in each 
r 

element of the residual matrix d„ and d < 12 <> 

f r 

(Results accurate to only 8 places). 

**If d is set to zero no residual matrix is punched. 

This can be an Illiac time saver if you feel the number 

of factors you are extracting are sufficient. However, if 

you wish to extract more factors at a later date you will 

need the residual matrix and the operator tape that follows 

it. To do this - follow "OPERATOR'S INSTRUCTIONS" above - 

but clear machine and read in operator's tape first. 

Machine will stop on 2^0SJ. At this point read in KSL 1.20 

master with a black switch up and follow steps 2, 3> and k 

making sure that your R-matrix in step k is the residual 

matrix of previous run. 

3. Communality tape - this tape consists of a series of v 

signed elements, scaled 10~ and terminated by an N. 

If ones are used they are entered as +1, any other number 

used has a leading zero, i.e., +0992* 

k* Correlation tape - the correlation tape is the lower 

triangular matrix output of K-8, K-9 or any other 

triangular matrix of correlations. This tape consists 

of v(v+l)/2 signed fractions scaled by 10°° . The diagonal 

entries must be present, but it is unimportant what 

these are, for they are replaced by the communalities .* 

TIME ESTIMATE 1 and l/2 minutes for input of master, parameter and < 

communalities. Time in seconds for input of R-matrix 
~ v(v+l)d n . 

loo 



r 1 

Time in seconds for output = — + ■ > '\ — ^-Id 

Calculation time/factor = 1 minute if v < 60 

2 minutes if N > 60 
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FF STOPS FF - location 026 (l) overfactoring) (2) underestimation of 

communalities | (3) R-matrix is non-Gramian| (k) machine 
errors 1 and (5) errors in data tape. Raise -white switch 
up and down. Variance table, R/ x and residual communalities 
will "be punched out. 
FF - location O^J drum transfer failure. 

NOTES * If you wish to obtain an R-matrix with .diagonals some 

value other than one - this may be done by setting f = 
and d f = in parameter, assigning a value to d . No factor 
extraction will be performed. The machine will read in the 
R-matrix, replace the diagonals by the communalities you 
read in on the communality tape, and punch out the same 
R-matrix with new diagonals . 

If cL. and d are set equal to in the parameter, and f is 
assigned some value (usually greater then the number of 
factors you feel are present), the machine will print out 
only the variance table. From this table particularly the 
eigenvalues - the number of factors to extract can be 
determined . 



KSL 1.21 
ESTIMATED CENTROID WITH ESTIMATED COMMUNALITIES 

PURPOSE This routine takes a correlation matrix in lower triangular 

form (K-2, K~8, K-9, and K-17 output) and breaks it down 
into a rectangular matrix of f factors (columns) and n 
variables (tests rows) where n is equal to the number of 
columns and rows in the input correlation matrix a The 
elements of the output matrix of 1.21 are called factor 
loadings » The highest valued element (+ or -) in each row 
of the correlation matrix is selected and replaces the 
diagonal element (usually +100) in the correlation matrix 
inputted. It will also give a residual correlation matrix 
after all factors are extracted. r .The; program will accept 
no more than 111 variables (tests). The number of factors 
extracted is limited to < n-1. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 1.21 (part l) 2^3F7 

2. Parameter 2^032 

5- Correlation tape 3^3F7 

4a. Master tape (resume at place removed 

for l) black switch up 2k^>F r J * 

kc Master tape, black switch down OF ** 

NOTE Master must be re-read and above instructions followed 

for each new set of factors extracted. 

* At this point only "f" factors signified in parameter 
will be extracted (no residuals). 

** By reading master tape into end, residual matrix is given, 

Step ka may be eliminated on operator's instruction sheet 

indicating k instructions . 

TAPE PREPARATION Master - Duplicate library master copy KSL 1.21 

Parameter - 004+ 

OOF OOfF 
OOF OOd.F 
OOF OOd^F 
OOF OOnF 
2450N 



where 

f = number of factors to be extracted (must be less than 

n below) 

d.= number of decimal 'digits in each correlation 
—l 

d= number of decimal digits in output 
—0 

n = number of variables (tests) or rows or columns in 

correlation matrix <» 
Correlation tape - In general this will be the triangular 
output from programs K=2, K-8, or K-9- It will be scaled 
10 3 signed (+ or -) and will have +100 in the diagonals, 
although any diagonal entry will be replaced by highest 
value in column » 
For examples +1000 

+0936 

+1000 or +1000+0936+1000-0032+0625+1000 

=0032 

+0625 

+1000 

For this example f could be 1 or 2| d. = h% d^ = 3 to 11 J 

n = 3- 

ESTIMATED TIME IN SECONDS 

Input of Master - 1 minute 

T + * A + 3v 2 d. + 13v 2 + 15000 
Input of data - 1 

100 

[f (12, 5v 2 + 17d + 150v +1000) +20,000] 
Extraction of factors - — . 

^. 1 . • = 1 E"(v)(v+l) /., u/,.] seconds 
Punching residuals ~ l s ' v - (d + 4)/50 J 

For example, if v = 10 3 f = 5 ^ and d. = k 

Then time 

_. . +1200+1300+15000 , [5(1250+68+1500+1000) +20,000] 

1 rain ° " 1606 + 1000 

+hhO r^ _ . , 
—=rr— = 2 minutes o 
50 



2 - 



FF STOPS 



NOTE 



FF location 148/ ^x. Error maybe due to 5 (l) R-matrix 
is non-Gramian| (2) machine error, and (3) errors in data 
tape» By raising white switch up and down, residuals 
will be printed o These may help in discovering source of 
error . 
The output from KSL lo21 is in the following forms 



CENTROID FACTORS 
ESTIMATED COMMONALITIES 
+0004.170 



+0..223) 
+0.315) 

+0000 c 960 



+001 +003 



Factor I 



(title) 



(absolute sum of elements 
in square correlation matrix 
inputted ) 



(absolute sum after 1 factor 
extracted) 

(reflections af terL:Ist 'factor ) 



RESIDUALS 
+0000.112 



(absolute sum in residuals) 



+0o010 ) 
+0i023 -0*032) 



residual correlations 



(operator tape to be used in further 
extraction of residual matrix at later 
date) 



00400+ ) 

80F 00001F ) 

OOF 00000F ) 

24999N ) 

As may be readily seen, to use the factors from KSL 1.21 

in KSL 1«91, KSL I08O or other routines using unrotated 

factors 1 the factors must be placed on a new tape with 

N's at end of each column » 

DO NOT PLAGE OUTPUT FROM KSL 1.21 BACK INTO ILLIAC WITHOUT 

MAKING NECESSARY DELECTIONS AND INSERTING AN N AFTER EACH 

FACTOR COLUMN o 
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KSL lo52 
ESTIMATION OF COMMUNALITIES BY METHOD OF MAXIMUM LIKELIHOOD 

PURPOSE This program -will calculate and print out a column of 

communalities that may be used in KSL 1.20 (Centroid Analysis 
with Fixed Communalities ) • It is limited to a correlation 
matrix (R-matrix) not larger than 8^ variables (Qk rows or 
columns in K-8 output);"': Through an iteration process 3 it 
converges to a relatively fixed level. The user must 
indicate how many factors to be used and how many iterations 
the program is to make. He also indicates at which point 
he wishes the iterations to be punched out. The output also 
includes a variance table^ giving number of factors ^ number 
of iterations and the amount of variance extracted by each 
iteration* 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 1.52 2402S 

2. Parameter 2^038 
3* Communality tape 2^0^N 
ho Correlation R-matrix 2402S 

(At this point j a new parameter _, communality tape^ and 

correlation R-matrix may be read in without re-reading 

master. Should the white switch be moved up and down, 

another iteration will be performed). 

TAPE PREPARATION Master - Duplicate library master copy KSL 1.52 

Parameter - k unsigned decimal numbers j a fifth hole 

character (i, : e», space,, carriage return or delay) separating 

each decimal character. The characters must be in this 

orders f space d space I space T space 

where 

f is the number of factors for which communalities are to be 
calculated . 

d is the number of decimal digits in output . 

I is the number of iterations before printout occurs • 
(This must never be larger than T below) 

T is the total number of iterations desired* 



COMMIMALITY TAPE 



CORRELATION R-MATRIX 



ESTIMATED TIME 



FF STOPS 
NOTE 



The number of positive signed values in this tape is equal 
to the number of variables (rows or columns) in the correla- 
tion R-matrix. This column of numbers must be ended with 
an N. The numbers may be l's (+100) if there is no idea 
of what the communalities may be. However 3 time will be 
saved on the machine if trial values less than one are used, 
These values are signed and scaled 10~ (i«,e<._, a value of 

.56 would be +056)0 Many times a trial communality tape 
will consists of the highest value of each row in the correla- 
tion matrix. 

In general 3 this will be the triangular output from YL-2, 
K-8 or K-9° It will be scaled 10~ , signed (+ or -) 
and will have +100 in the diagonals » 

fTF[l3v 2 + 36v + 55] + T[2v 2 + l8v + k6o] + (T-l)(l75v) 
+1500*) seconds x 10~ 5 . 

For example j a 10 factor | 30 variable; 20T| 101 problem 
would be £200 [11700 + 1080 + 55] + 20 [1800 + 5^0 + 460] 
+52500 + I5OO J 10~ 5 or = ^5 minutes 
(See parameter description for codes v = variables) 
There are no FF stops listed for this program . 
Output is in this form 
+0239 
+0327 
+0*02 
N 



+0241 
+0326 
+0^33 
N 

VARIANCE TABLE 
The particular iteration to be used in KSL 1.20 can be 
duplicated from output tape. Use only the particular < 
column desired plus the N following it as a communality 
tape for KSL 1.20. 
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PURPOSE 



KSL 1.53 
EXTRACT ELEMENT WITH LARGEST ABSOLUTE VALUE FROM A MATRIX 

This routine takes the lower half of a triangular matrix 
(K-2, K-8, K-9 output) and extracts the largest absolute 



value element, printing it out in columnar form, scaled 

and signed with +'s. An N terminates the column - the 



10 



output is thus available for use as a eommunality tape 
in KSL 1.20. The size of the matrix (N) must be < 900. 



OPERATOR'S INSTRUCTIONS 



1. Master tape 
2 • Parameters 
3. Matrix tape 



Stops 
3^00S 
24027 
2400S 



PREPARATION OF TAPES 



TIME ESTIMATE 
FF STOPS 



A new parameter and matrix tape may be read in without 
rereading master. 

1. Master - duplicate library copy KSL 1,53 

2. Parameter - d Q space n space where d = number of decima 
places in output < 9 and n = size of matrix 

3. Matrix tape - this consists of n(n+l)/2 signed fractions 
in triangular form as given in K-2, K-8, K-9. 

Time in seconds = v (.002d + .006) + v(.019d n + .053 ) 







FF000 location 0^3 - sum check failure on master. 
FF001 indicates that d Q in parameter > 10. A new parameter 
; can be read in with white switch up and down. 
FF002 indicates a sign digit is out of place. White switch 
up and down will cause matrix to continue to be read. 



KSL 1,71 

OBTAIN CORRELATION COSINE MATRIX AMONG PRIMARY FACTORS 
AND FACTOR PATTERN STRUCTURE 

PURPOSE Given a cosine matrix between reference vectors (factors) 

and an oblique factor matrix of reference vectors , this 
program will provide the cosine matrix (symmetric) among 
primary factors and a factor pattern matrix of the same 
number of factors and variables as the inputted factor 
matrix of reference vectors. The size of the cosine 
matrix fed in and received as output must be no larger 
than 22 « The number of columns (factors ) in the reference 
vector matrix must be equal to the number of columns in 
the cosine matrix. Generally ^ the V ^ matrix (primary factor 
pattern) obtained in this program is used for KSL 1.72. 
The two matrices obtained from this program should be 
identical (except for rounding error) to the corresponding 
matrices obtained from KSL 1.91« 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL I.71 2^3F7 

2. Parameter 2^035 

3. C tape 3^0F7 

k-. V tape OF 

rn * 

NOTE Master must be re-read for each problem. 

TAPE PREPARATION Master - duplicate library copy KSL I.71 

Parameter - 0010+ 

OOF OOf F 

OOF OOd.F 

1 

OOF OOvF 
2^53N 
where 

f = number of factors or columns in reference vector matrix 
d. = total number of decimal places in one variable of 

reference vector structure (< 12) (all variables contain 
same number of digits) 



v = number of rows or variables in reference vector matrix 
C tape - punch a square symmetric matrix. Each element is 
signed and each element should contain no more than 11 digits. 
Terminate each row with an N. Scale each element 10" . 
V tape - Punch by rows - no terminating symbols. Each 
variable of element should be signed and contain no more than 
11 digits . 
NOTE In the V tape any element of +1.00 must be entered as a +099 < 

ESTIMATED TIME C f : order of C f time (in seconds) 



5x5 


ko 


10x10 


70 


15x15 


90 


20x20 


120 


22x22 


1^-0 


NO. DIGITS x v x n 


^o 



V time = 50 — 
where 

v = number of variables or elements in row of V 

— rn 

n = number of rows in V 

— " rn 

FF STOPS FF hangup in location 02J indicates a drum failure. Re- 

read master, try again. 
Hangup on ON3276NN3277 in location 139 / l6 ) indicates that 

matrix C has no inverse, 
r 



KSL 1.72 
ESTIMATION OF OBLIQUE FACTOR SCORES 



PURPOSE 



Given four matrices: F^ 



factor matrix: Vp 
' -fp 



arbitrary unrotated orthogonal 
- primary factor pattern obtained from 



KSL 1.71: -A- - matrix that transforms F. to V : and 
P * — r 2 -0 -m* 

£ - a single column of 1-h (see preparation of tapes 

for description) this routine gives the three matrices: 

—A- - a transformation matrix that takes F~ into V- : 
— f -0 -fp 5 

V_ - the primary factor structure matrix and; B - an 
n x k matrix of regression -coefficients. The B matrix 
can be multiplied by the Standardized Raw Score Matrix to 
obtain oblique factor scores 



The matrices F». V- . 
-0* — fp* 

■V- must have no more than 22 columns (F) or 111 rows (v) 



OPERATOR'S INSTRUCTIONS 



Stops 
2404F 
3*028 
3^051 
k. _/u tape 2^154 

5. F n tape 2404F 

(A new problem may now be run by inserting a new l-h^ 
tape and proceeding through steps 3.> ^j and 5)« 
PREPARATION OF TAPES Master - duplicate library copy KSL 1.72 



1. Master KSL 1.72 

2. 1-h tape (residual communalities ) 

3. V fp tape 
/u tape 



1-h 



This is a column of unsealed fractions. An N 



must be punched at the end of this column. None of these 

2 

fractions can be smaller than +.01. h is obtained 

by summing jthe squares of each element in a row of F~ . 



EXAMPLE 



h 1-h tape 

1.13 +01 (see above) 
.90 +10 
N 
output tape from KSL 1.71. 
It may be copied directly from 1.71 output arid used. 



+.90 +.30 -.20 
+.20 + 99 * 30 

+.4o + 50 + 70 



V- tape: This is the V 
-f p * — f p 



It is punched row by row. Each row has an N at the end. 
Each element is scaled 10 



JL ; This transform matrix is punched by columns | scaled 

-1 
10 and has an N at end of each column. If the transform 

is used from KSL 1.9&* it must be remembered that it is 

output from KSL I.96 by rows. Therefore, it must be 

transposed or punched by columns. 

F^ tape: This is punched row by row. Each row has an N 

at endo Each element should be punched as a signed s 

unsealed fraction of no more than 11 digits . (This will 

generally be a matrix of factors obtained from either 

KSL 1.20 or KSL 1.21. The results from these two programs 

must be transposed (see KSL 5»30) and scaled 10 before 

using on KSL 1.72). 

TIME EST1MA.TE None available 

FF STOPS FF location 02J - drum failure - re-read master 

FF location 1^5 - wrong number variables in a row 
0N3276M3277 location 139 of V - -A_ is singular 

•*" JET ■*- 

NOTE The rows of -^ = columns of ^a- (-£- must be square ) 

The rows of F^, = rows of Y\, = the number of elements in 
2 -0 -fp 

1-h tape 

The columns of F A = columns of V„ = rows of-*" . 
-0 -f p — -r 
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KSL lc80 
VARIMAX ROTATION 

PURPOSE This routine takes a set of unrotated factors (a f x v) 

matrix (usually obtained from KSL 1.20 or KSL 1.2l) and 
rotates them maintaining an orthogonal reference system 
according to the Varimax criterion. Each output contains 5 
kinds of data in this orders (l) the communalities (h 
see KSL 1.72) of the v variables, scaled at 10 " 2 ; (2) a 
set of criterion values of 12 digits each, scaled 10 . 
These values are used by the machine to determine when con- 
vergence is reached! (3) the variance contributed by the f 
varimax factors, scaled at 10 " | (k) the parameters v, 
f , and d~ (these will be identical to parameter inputted 
(see below) )j (5) the £ variamx factors scaled 10 . 
Each column or f varimax factor is followed by an N. 
The number of variables must be < 127* The number of 
factors must satisfy this inequality : (N+l)(f+l) < 102^0. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 1.80 2400K 

2 . Parameter 240^3 

3. Data (black switch down) OF 

(To read in new parameter and data, center black switch 
and white switch up and domi. Repeat 2 and 3 above). 

(if problem must be removed before OF is reached, center 
black switch! when 3^0SF occurs white switch up and down. 
Output will ensue and parts (k) and (5) of the output tape 

(see PURPOSE above) may be then run at a later date to 
continue rotation problem). 
PREPARATION OF TAPE Master - duplicate library copy KSL 1.80 

Parameter - v space f space d_ space in that order, 

where 

v = number of variables 

f = number of factors, and 

d~ is the number of digits desired in output '(since output 

-1 
is scaled 10 , a d of 3 will contain only 2 significant 

figures ) . 



Data - This is an unrotated factor matrix punclied column 
by column. It is sealed 10° s has a sign in front of 
each element and has an N at end of each column. Output 
from KSL 1,11, KSL 1.20, and KSL 1.30 may be used as 
data without changes » KSL 1.21 output may be used if 
correct changes are made (see KSL 1.21). 
TIME ESTIMATE 200 + .OO^vf 3 seconds is a conservative estimate 
FF STOPS None given other than library subroutines. 



KSL 1,90 
OBLIMAX ROTATION OF FACTORS 

PURPOSE This routine takes an f (factor) x v (variable) unrotated 

matrix and rotates it to an oblique set of factors . See 
the library printup of this routine (pages 4-7) to obtain 
an idea of the mathematical procedure used^ and its limita- 
tions o According to options specified on parameter tape 
(see below) the user may obtain varying combinations of 
results o The limits to this program are given under 
NOTES . With appropriate manipulations to the output , it 
can be used in KSL 1.95, KSL 1*96 or KSL 1.72* See those 
routines to determine what manipulations are needed. 

OPERATOR'S INSTRUCTIONS Stops 

A. 

1, Master KSL 1.90 2407J 

2o Parameter tape 240F0 

3„ Factor tape^ run with black switch down OF 

If black switch is placed in center _, machine will stop 

at 3^-092 after each iteration . A white switch up and 

down at this ^>k stop will cause print out of results 

obtained at that point , This allows interruption of 

problem before it is completed. To read another parameter 

tape and factor tape once OF is reached 3 white switch up 

and down 

Bo Subsequent run with a transform (either from interrupted 

run as in A or if transform is available at beginning of 

run). 

Stops 

lo Master KSL I08O 2407J 

2o Parameter tape 2^0F0 

3o Transform tape (white switch up and 

down) 3^0L8 

4o Factor tape,, run with black switch down OF 

PREPARATION OF TAPES Master - duplicate library copy KSL lo90 

Parameter - This tape consists of four unsigned integers ^ 

each separated by a fifth=hole character (carriage : 



return-space-delay). They are d space v space f space and 
x; where d„ is the number of digits in. output (results are 
scaled 10 - ask for one more digit than you wish significance 
in); v is the number of rows or variables in one factor; f 
is the number of factors or columns in the factor tape, 
and x is a directive telling what options to print out. 
(See below). 
OPTIONS A directive may run from to 15 . 

Directive Results obtained 

criterion values only 

1 T and V (by columns - used for 

subsequent runs) 

2 T and V by rows 

r J 

3 T and V by both rows and columns 

k C , n, V^D" 1 and D 

t 3 f 3 i 3 

5 T and V by columns C , C f , V 

D"* 1 and D 

6 T and V by rows ; ■ C . CL , V- 

D and D 

7 T and V by rows and columns ; 

C t 3 C f 3 V f > D_1 and D 
Increase thru 7 ^y 8 (i.e., 8 thru 15 ) and the bblimax 

rotation process will be suppressed, [if you have a T 

and V rotated and wish to get C^ and V a directive of 

ik will do this without the rotation process being performed 

in the machine] . 

In the above directive options the following symbols are 

useds 

T is the transformation matrix (scaled 10 ) used to transform 

the factor matrix to a 

V which is the rotated factor matrix (scaled 10 ) 

C is the reference vector correlations . It is obtained 

r -1 

by finding T'T, is symmetric, scaled 10 and is punched 

in triangular form by rows 

C_p is the primary factor correlation matrix punched in 

-1 
triangular form by rows and scaled 10 

_ 2 - 



V- is the primary factor pattern punched by rows, N at 

-1 
the end of each row, J at the end of matrix, scaled 10 

D is a diagonal matrix of the reciprocals of the square 
roots of the diagonal elements of C " . D is punched 
as a column is unsealed and has UJ at end. 

D is the inverse of D. Decimal points indicate the 
scaline. It has an NJ at end of column. 

THE TRANSFORM TAPE 

(For a run as described under operator's instructions, B, 

above). Punch by columns with an N at end of each column. 

They must be scaled 10 and have signs. 

THE FACTOR TAPE 

A set of factors (obtained from KSL 1.20, 1.21 or other 

factor program) punched by columns with an .N at the end of 

each column, signed fractions and scaled 10 
TIME ESTIMATE 1. Read master - 2 minutes 

2. Read data tape - v f (.009 + .OO^d.) seconds 

3. Approximate calculation time/iteration « 
f 2 (.230 * .007f + .010v) seconds 

(The number of iterations varies for each problem. The 
larger the problem - the more iterations ) . 
k. Printing time when x = 5 or 6 

[5f + *35^f 2 + .012f 5 + .033d Q [f 2 + vf]] seconds. 
For example, a v of 50, f of 10 and d. and d n of 5 might 
take 10 iterations. Time would be approximately 2 minutes 
+[500 (.009 + .016) + 1000(o230 + .07 + .5")] 
+ [50 + 35A + 12 + 19-6] = 18 minutes 
(Since the number of iterations are undeterminable for 
large problems, it is suggested that the user run a problem 
for a block of time, after which the Illiac operator is 
instructed to remove problem as described above in Operator's 
Instruction. The problem may then be finished in a subsequent 
run or two. To do this, both T and V are required by 
columns for subsequent machine runs. The directive x (see 
PREPARATION OF TAPES - parameter) should be an odd number). 



FF STOPS FFCOA - sum check failure - reread master 

FFOOS - C matrix is singular (has no increase) 
FF015 - limit of factors exceeded (see NOTES below) 

NOTES If the directive x (see parameter tape) is k, 5 ,.-6, 1 9 

12, 13, lh or 15 the limit on the number of factors, f, 
must be less than 21 , For any other x directive the limits 
are d Q < 12; v < 128, f < HO. 



KSL 1.94 (Procrustes) 
SOLUTION OF MATRIX EQUATION At=B WHERE A AND B ARE KNOWN 

PURPOSE Given two matrices A and B, this program will give you a 

Transformation matrix (if one exists ) which will transform 
matrix A into matrix B| it also punches out At. The columns 
of A and B < 22. The rows of A or B < 136. The number of 
columns of A need not equal the number of columns of B 
but the number of rows in A and B must be equal. The sum 
of squares of columns of t are equal to one (the columns 
are normalized). Because of this the value of At will not be 
exactly that of B. 



OPERATOR'S INSTRUCTIONS 



1. Master KSL I.94 

2. Data A 

3. Data B 

4 o Parameter 



Stops 
24046 
3409N 
340S1 
24046 



PREPARATION OF TAPES 



TIME ESTIMATE 
FF STOPS 



A new problem may be begun by reading in a new data A and 

following steps 3 and 4. 

lo Master - duplicate library copy KSL 1„94 

2 Data A - this is punched by columns with an N terminating 
each column. The last N is followed by a J. Each 
element begins with a sign^ has less than 12 digits and 
is scaled 10 

3. Data B - punch exactly as Data A 

4. Parameter - d n space where d^ is one character from 2 
thru S to indicate the -1 number of decimal places wished 
in output. 

No time estimate given. 

FF order indicates wrong number of elements in a column 

(as compared to the first column). 



KSL 1=95 
PLOTTING FACTOR PAIRS 



PURPOSE 



This routine films the plots of all factor pairs of a set 
(no. frames = * ' ). It would "be used if you only wish 
one set of plots 9 i^e.^ you do not plan on rotating further, 
If you do plan to do so 5 it is advisable to use KSL \>S6 3 
(one exception •= this routine would he used should the 
number of factors you wish to film exceed the capacity of 
KSL lc96)o 



OPERATOR ' S INSTRUCTIONS 



PREPARATION OF TAPES 



TIME ESTIMATE 



FF STOPS 



Master KSL 1„95 



Parameter 

Set output to scope 

Data (black switch down) OF 

Clear machine - advance camera one frame 



Stops 

2^3F7 
2^267 



k ID 

1. Master - duplicate library copy KSL 1.95 

2 o Parameter 00279-h 

OOF OOvF 
OOF 00(f-l)F 
2^267N 
where v = number of variables 

f = number of factors 
3o Data - this tape consists of the data you wish plots 
on« It is scaled 10 has signs in front of each 
element,, an N after each column "and a J after last No 
k ID tape - this tape consists of any identification you 
wish placed on film to help you identify it. The 
first character is a 9° Your identification follows. 
A one-hole delay follows the identification. Do not 
put a one-hole delay in the identification . 
A very rough estimate can be obtained by this formula . 
Time = 2f (f+l) seconds 
Only in library subroutines . 



KSL I..96 
ROTOPLOT ROUTINE 

PURPOSE Given an unrotated (Fxv) matrix (obtained usually from 

KSL 1»20 or KSL 1.2l)| a transformation matrix (fxf ) 

and a shift matrix (fxf) this routine will multiply the 

8 matrix times the T matrix, normalizing columns in the 

resultant new T. matrix - multiply this new T. matrix 
1 i 

times the unrotated matrix to obtain an oblique V matrix o 

rn 

It outputs the T. matrix, the V matrix, a hyperplane 

count (see below) and the upper half of the cosine of 

angles matrix between the reference vectors (f)» The 

routine is limited to m variables and 22 factors o The 

output of the X and V matrices is unsealed and limited 
.f rn 

to 2 decimal places » The routine then places on film 
scatter plots between all factors; (l with 2, 1 with 3* 
1 with .k 9 etc. 5 2 with 3, ...) until A / frames have 
been filmed » This program is used when you wish to 
visually rotate the output to arrive at simple structure o 
OPERATOR'S INSTRUCTIONS Stops 

lo Parameter (preceded by DOl) 243F7 on one 

.2. Jfester tape KSL lo96 2kO?S ^J^ 

3o Data .1 (S. matrix) 2^0NK tion of 

ko Data 2 (\. matrix) 240N8 tapeS ^ 

.1 

5» Data 3 (V Q matrix) 2^03N 

When Illiac stops on 2^3^ at end of Data 3 set output 

to scope o 

60 Black switch down OF 

When OF is reached, advance one frame, clear Illiac and 

read in ID tapeo 

NOTE On operator's instruction sheet to indicate the number of 

f ( f +1 } 
frames required use formula v, +1 where f = number of 

factors o 

PREPARATION OF TAPES Master tape - the master tape can be duplicated from library 

copy» It is preceded either on same tape or by extra 

smaller tape to be read in first « This extra tape contains 

a copy of the DOI (in library) and the parameter (see below) „ 



Parameter tape 0010+ 

00F OOrF 
OOF OOf F 
OOF OOvF 

where r = number of rows in transformation tape (X. tape) 

f = number of factors in V^ tape 
OOF OOvF = number of variables in v' tape 
Data 1 (S. matrix) 

This is a shift matrix (to be multiplied times the previous 
transformation (x) matrix to give a new X matrix). It is 
punched by rows,, with +1000 in the diagonals and the other 
elements signed and scaled HT . Each row is terminated 
by an N. The final N is followed by a J. On the first 
run of this program,, the shift matrix is usually an identity 
matrix,, i.e., +1000' s in the diagonals and 10000' s for 
every other element. A + can be punched for a +0000. 
Data 2 (X. matrix) 

This matrix is punched by rows <> Each row is terminated 
by an N. The last N is terminated by a J. Elements are 
signed and unsealed. The first run on this program re- 
quires a X matrix punched by hand. It is usually obtained 
from any KSL routine that gives a X matrix. On subsequent 
runs with the same V n , each new X matrix can be duplicated 
directly from the output of the preceding run*. (Do not 
duplicate the heading on the output tape ) . 
Data 3 (V n matrix) 

This is generally the output from KSL 1.20 or KSL 1.21. 
Since the Data 3 is unsealed and punched by rows 5 before 

using output from KSL 1.20 and 1.21, it must be run on 

+1 
KSL 5° 30, transposed and scaled 10 . This will automati- 
cally complete the other requirements , an N at end of each 
row 5 a J after final N. 

ID tape - this tape places on film your name and any 
identification you wish. It is constructed by punching 
a 9 and then what you wish to appear as identification. 



A one-hole delay comes at the end , and xmLy the end of 
your identification o 
TIME ESTIMATE Time = ~ (f + g + 9) minutes 

where f = number of factors 

v = number of variables 
FF STOPS FF indicates wrong number of variables in a row. 
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KSL lo97 
MAXPLANE (An Oblique Factor Rotation Program) 

PURPOSE This routine takes in a transformation matrix and a signed 

V n centroid matrix obtained from KSL 1.20 or 1.21 performs 
a set of oblique rotations in an attempt to arrive at 
simple structure through use of criteria introduced by the 
user and included within the program*. When these criteria 
are met, the machine punches out the rotated matrix, the 
new transformation matrix and the correlation matrix of 
angles between the factors. The rotated factors are then 
plotted on film (l vs 2, 1 vs 3, . .«, 2 vs 3* 2 vs k, „„,,, 
etc . ) until p — - factor pairs have been plotted « The 
number of factors < 20 | the number of variables < lil» 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL lc97 20— 

black switch up 2h location 001 

■2, Data tape (black switch down) - read in order 

A new set of data may be read in without rereading master 

by inserting the new data in the reader - while the previous 

problem is calculating • If during a run, you wish to 

interrupt the run, center black switch . The machine will 

stop on a 20 at end of a cycle. If the white switch is 

moved up and down .- the first three elements -of ."the output 

data (see below - Output Data Form) will be punched and 

the computer will stop on a 2k order . If the black 

switch is placed down .- computation will begin again . If 

the white switch is moved up and down .— a complete data 

output will ensue., 

PREPARATION OF TAPES 1. Master - duplicate library copy KSL 1*97 

2 o Data tape - the data tape is rather complicated o An 

error on the data will cause the machine to stop on 

an FF stop (see FF stops). The data tape consists of 

the following - in order given 2 

d Q s number of digits after decimal point in output (+, a 
number 1-9* and an N) 



t, s greatest hyperplane thickness 

tpS middle hyperplane thickness 

t s least hyperplane thickness 

(+, 0, decimal pointy thickness as an n place decimal 
fraction < 10 for each of the three thicknesses . The 
third thickness is followed by an N). 

X~i an identity matrix (+ ones in diagonal, +0's elsewhere) 
if a centroid V matrix from KSL lo20 or 1,21 is used, 
or the output matrix from KSL 1.10 (max. 20x20 ) (sign 
(+0»), (i,e.j 10"-'-) element value as an n place decimal 
fraction < 11 for each element . Elements are punched by 
column - each column terminated by a J. Last J followed 
by L) 

V Q s V matrix output of KSL 1.20 or KSL 1.21 always (max. 
111x20 )o (Sign, (i.e»j 10" 1 ) element value as an 
n place decimal fraction < 11 for each element in matrix). 
Elements are -punched by columns . Each column terminated 
by a J. The last J is followed by an L. 

Example of data tape? 3 factors | k variables 

d Q = +5N 

t x = +0.10 

t 2 = +0.075 

t^ = +0.0125N 
3 

+1000 +0000 +ooooj +0900 +0500 +0800J 

X = +0000 +1000 +0000 J or +0001 +0700 +0600J 
+0000 +0000 +1000JL +0300 +0510 +0000JL 

+0923 -0621 +0321 -021U 
-0326 -0921 +0032 -0022J 
+0^32 -0222 =0662 -0331JL 
FORM OF OUTPUT; 1. The first integer refers to the number of rotations 

performed on data . 
2o A 3-column table giving the factor number, number of 
variables in that factor that lie within a + or - 
value equivalent to least hyperplane thickness (in 
example above - those variables falling between + o0125)e 
3. A correlation matrix of cosines between factors by 
columns . 



4. The parameters d„^ t, s t , t * (to be used for future 
runs with X data output should problem need to be 
interrupted and removed from machine - see Operator's 
Instructions ) <> 

5. X t a new transformation matrix formed by the routine., 

y n 

(Same form as \~ except with d n decimals). 

6. V : this rotated matrix has same form as V~ (with d~ 

n 

decimal places) and is obtained by multiplying V\A » 

7» Scope plots for each different pair of factors , plotted 

on orthogonal axes . Each axis is labeled with its V 

n 

matrix column number. The upper right hand corner con- 
tains the associated C matrix entry. 
NOTE * 5 carriage returns and 15 two-hole delays are: printed 

between the last element of the C matrix and the d~ on output 
tape to enable the user to find and read required portion 
of output data tape back in to Illiac should the occasion 
arise. 

TIME ESTIMATE 1. Read master tape - 2 l/2 minutes 

o 
2. Read input data in seconds .= .07 + .004 f (d.) 

+ .604 fvd. 

l 

3« Calculation time - variables f < 5 J v < 50 run less 
than an hour. Larger problems should be run for a 
fixed time - removed as instructed in Operator's 
Instructions ~ and rerun from this point at a later 
time. If rerun ^ be sure to use outputted X with V*. - not 
with outputted V . 

ko Punch and plot output in seconds i 

£ (4+d ) 

)0 — f (f+v) + 3(f-l) (f ) 



3 



FF STOPS FF1; Improper terminating symbol on d Q 

FF2; Improper terminating symbol on t, t 9 t 

FF3l Incorrect number of t's input 

FF^-s Improper terminating symbol on X, n 

FF5s Incorrect number of elements as a column of X~ 

FF6s Incorrect number of columns input for X, n 



FF75 Improper terminating symbol on V Q 

FF8; Incorrect number of elements in a column of V Q 

FF9: ^Incorrect number of columns in V Q 

FF10: Drum error 



- k - 



KSL 1.98 
ROTATION MASTER 

PURPOSE Given an unrotated (Fxv) matrix (obtained usually from '.: 

KSL 1.20 or KSL 1.21); a transformation matrix (fxf) 
and a shift matrix (fxf) this routine will multiply the 
S matrix times the T matrix, normalizing columns in the 
resultant new T. matrix, multiply this new T. matrix 
times the unrotated matrix to obtain an oblique V matrix, 
It outputs the T. matrix, the V matrix, a hyperplane 
count (see below) and the upper half of the cosine of 
angles matrix between the reference vectors (f). The 
routine is limited to 128 variables and 52 factors. The 
output of the X aM V matrices is unsealed and limited 
to 2 decimal places. The routine then places on film 
scatter plots between all factors; (l with 2, 1 with 3, 
1 with k, etc.; 2 with 3, ...) until l n ^ n ~ ^ -frames have 
been filmed. This program is used when you wish to 
visually rotate the output to arrive at simple structure. 

OPERATOR * S DESTRUCTIONS Stops 

1. Clear start parameter tape 2^05J 

2. Master KSL I.98 black switch up 2*4-188 

3. Data 1 black switch up 2^163 
k. Data 2 black switch up 2*1-172 

5 . Data 3 black switch up 0F000 

6. Clear start label tape 
PREPARATION OF TAPES Master - duplicate library copy KSL 1 . 98 

Parameter tape - 009+ 

OOF OOeF 
OOF OOpF 
OOF OOdF 
OOF OOrF 
OOF OOfF 
OOF OOvF 
2*J-999N 



where 

e* =? number of elements on a line of the print-out 1 < e < 70/(p+l) 

p = number of significant places in the output 1 < p < 12 

d = dimensions = number of rows of X matrix 2 < d < 52 

r = number of rows of S matrix 2 < r < 52 

f = factors = number of columns of S matrix 2 < f < 52 

v = variables = number of rows of V~, V matrices 2 < v < 128 

0' n — — 

Data tape - data tapes are prepared precisely as for KSL 1.^6. 

Output - output is in the same order as with KSL I.96. 

However, there are slight differences: 

1. With large matrices, it is impossible to print an 

entire row of the matrix on one line of the teletype 
paper . This program prints e elements on a line, then 
executes a carriage return and line feed, (e is a 
program parameter). At the end of a matrix row, a 
carriage return and two line feeds are executed. 
(Exception: spacing of the C matrix is twice this e.g,, 
two line feeds between lines, and four line feeds between 
rows ) . 

2.-:.. Output is not \ ■'limited to two significant figures. Instead, 
the number of significant figures is determined by a 
program parameter, p. If p=5^ l/l6 would be punched as 
+06250j if p=2, l/l6 would be punched as +06. 

3» All matrices are punched by rows. The C matrix is 

punched as an upper triangular matrix, with the diagonal 
entries (which are all one) deleted. Due to restrictions 
placed on the format by the length of a teletype line, 
the C matrix will not appear as a triangle. Nevertheless, 
it is ai upper triangular matrix. 

h. On the oscilloscope output, a "+" in the extreme upper 

right corner is to be interpreted as a zero, and a one 

to carry, i.e., : 

+57+ means +58O +09+ means +100 

7-09+ means -100 
* e should be < 17 if p = 3 



FF STOPS 
TIME ESTIMATE 

NUMBER OF FRAMES 



A time estimate is given in terms of the program parameters 

f, v } and p. The number of frames is given in terms of the 

program parameter f . 

WARNING : Any attempt to use parameters &_, r, f, or v larger 

than the maxima sho"wn may cause Illiac to hang up. Even if 

Illiac does not hang up^ the results are guaranteed to be 

incorrect . 

FFOOO^? Both denote drum transfer failures. FFOOO is a 
FFOlOj program stop, FF01G is the Yl subroutine stop. 

220 +l) + fo = k minUteS 
where f, v, and p are defined under Parameter Tape Preparation. 
1/2 f(f-l) + 6 
where f is the number of factors. 
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KSL 2,01 
K-8 REVISED TO ACCEPT SINGLE, UNSIGNED DIGITS 

PURPOSE This routine will give you the same results as obtainable 

in Illiac library routine K-8 (Pearson Product Moment 
Correlations ) o The number of variables is limited to ikk, 
The sample size or number of observations is unlimited e 
The output contains columnar correlations followed by 2 
columns, the first means, 2nd standard deviations. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 2*01 3^01K 

2, Parameter 2"+0S9 
.3. Data 2^01K 

A new parameter and data tape may now be read into Illiac , 

PREPARATION OF TAPES 1. Master - duplicate library copy KSL 2,01 

- A 

2» Parameter - -cLF d-L dLJ ssS -, _ 

1=0 



,-!■ 



where d n F = decimal output in correlation (results scaled 10 ) 
d n L = decimal output in covariances (results scaled 10*" ) 
d n J = decimal output in means and standard deviations 
ssS = number of rows of observations! sample size 

= triangular output (use this if you wish to put 

results back in Illiac. Do not get covariances 
if you are reading results back into Illiac, 
i e«, set d~L = 0L« 

1 = square output (symmetric) with N's.at end of 

each row or column . Use only if you wish to go 

on further with output on Illiac. 

3<. Data tape - this consists of a set of single unsigned 

digits (from 0-9) .with sin N terminating each row of 

variables . 

ESTIMATED TIME This may be obtained from the time chart included with 

library writeup of K-8, 

FF STOPS FF015 = wrong number of variables in a row 

FF0l6 = variance is negative 
FF017 = sum check fails on master 
FF023 = arithmetic error, possible from scaling. 



KSL 2.03 
MEANS, STANDARD DEVIATIONS, THIRD AND FOURTH MOMENTS ABOUT THE MEANS 

PURPOSE For each of a set of v variables , this routine will calculate 

the mean, standard deviation, and the third and fourth 
moments about the mean.. The program reads the data in the 
form of signed .multi-digit elements or in unsigned single 
digits 1 to 9= The number of variables is limited to < 1^-5° 
If means and standard deviations are required (for example, 
to be used in KSL J+.10) the program can be used with as 
much as a 9&fo saving in Illiac time over K-8 or K-0_7<> 

OPERATOR'S INSTRUCTIONS Stops 

1. -Master KSL 2.03 3^084 

2. .Parameter 2^0*fN 

3. Data tape 2^06^- 

A new parameter and data tape may be read in without re- 
reading master. 

PREPARATION OF TAPES 1. Master - duplicate library copy KSL 2„03 

2o Parameter tape - this consists of rNd JnFpL 

where 

r = if data consists of unsigned digits 

i\ ■■-■- 1 if data consists of signed digits 

d n = number of decimal places in M - output. The third 

and fourth moments punched to 10 places 
n = number of rows = i.e., sample size 
p = if no third and fourth moments are desired 
p = 1 if third and fourth are desired 
3. Data tape .- the v variables or tests are punched by 
rows and consists of n rows. Each row is terminated 
by an No* If r = 1 (see Parameter) decimal digits for 
each element is < 12. 

FORM OF OUTPUT Means and standard deviations are printed out in parallel 

columns terminated by an N (ready for use in KBL k .10).** 
If moments are calculated, they follow in like parallel 
manner, terminated by an N. 



TIME ESTIMATE Master tape readin - 30 seconds 

Data readin - .004 (d. + l) vn seconds 
Calculation - .003 vn seconds 

Punch - ,035v (d Q + l) for means and standard deviations 
•kv for third and fourth moments 
FF STOPS FF003 sum check failure 

FF000 location 0S1 wrong number of variables as compared to 

row 1 
FF002 location 0L5 variance is negative due to rounding 
errors and outside tolerance limit (l0~ 10 ) white 
switch up and down to bring about punching. 
NOTES .* If an N in data tape is replaced by an F, the machine 

will stop on a 3^ order allowing a new data tape to be 
read in with a black switch. This allows the user to split 

a large data tape in to smaller, more manageable data tapes. 

P 

** If the scaling on variable j is 10" j thenjin the results 

the means and standard deviations will also be scaled by 
10 j. The third moment will be scaled by .ICT^ j and the 



" ;J;P 
fourth moment will be scaled by 10 j 
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KSL 2.40 
PHI'S OR COVARIANCES FOR DICHOTQMOUS DATA 



PURPOSE This program reads in a data matrix consisting of responses 

coded 1 and and punches out a matrix of correlations 
or covariances (< .25 for covariances ) <> The number of re- 
sponses is < 157 • The number of persons or sample size 

<2535- 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL 2.40 ^ko66 

2. Data tape 2^1K4 

3. Parameter tape 2^-lK^ 

To read additional parameters to operate or save data at 
the end of step 3j raise hlack.'' -.switch*. To begin a new 
problem with different data at the end of step 3, raise the 
white switch up and down. 
PREPARATION OF TAPES 1. Master - duplicate library copy KSL 2.40 

2. Data tape - this consists of a series of ones and zeroes 
(l indicates a positive response , a negative response) 
These digits are unsigned. An N follows each person's 
response. The final N is followed by an L. 

3. The parameter tape is as follows s 

d Q space x space where d is the decimal output desired 
and x is a directive indicative of desired output. The 
meanings of the directive are as follows: 



X 


Coefficient 


Form 
Triangular 


Printed 
By rows 


Matrix 
Complete 







Phi's 




1 


Phi's 


Triangular 


By 


rows 


Submatrix 


(See below) 


2 


Phi's 


Square 


By 


rows 


Complete 




3 


Phi ' s 


Rectangular 


By 


rows 


Submatrix 




k 


Covariances 


Triangular 


By 


rows 


Complete 




5 


Covariances 


Triangular 


By 


rows 


Submatrix 




6 


Covariances 


Square 


By 


rows 


Complete 




7 


Covariances 


Rectangular 


By 


rows 


Submatrix 




8 


Phi's 


Triangular 


By 


columns 


Complete 





X Coefficient 

9 Phi ' s 

10 Phi ' s 

11 Phi ' s 

12 Covariances 

13 Covariances 
Ik Covariances 
15 Covariances 



FOEM OF OUTPUT 



TIME ESTIMATE 



Form 
Triangular 
Square 
Rectangular 
Triangular 
Triangular 
Square 
Rectangular 



Printed 
By columns 
By columns 
By columns 
By columns 
By columns 
By columns 
By columns 



Matrix 

Submatrix 

Complete 

Submatrix 

Complete 

Submatrix 

Complete 

Submatrix 



Note; If a complete matrix is desired,- the directive must 
be an even number - and the parameter is as given above 
(d n space x space )<» However, if an incomplete (submatrix) 
matrix is desired, the directive must be followed by four 
more numbers | i space g space r space s space where i and j 
are the subscripts for the first coefficient to be printed 
and r and s are the subscripts for the final coefficient 
to be printed. For example, if the Phi's calculated could 
form a matrix as such 

a ll a 12 a i3 S l^ 

a 21 a 22 a 23 a 2k 

a 31 a 32 a 33 a 3^ 

%± % 2 a k3 a ^ 
and you are interested only in the coefficients forming 



then your parameter might be d space 



this submatrix c. a 2 ^ a . 

a 33 & 3^ 
s (odd) space 2 space 3 space 3 space k space . 

If matrix "priced -is : ' triangular ; > - no terminating symbols. 

If the punched matrix is square or rectangular, an N will 

be punched at the end of each row or column . Phi coefficients 

are scaled 10* . Following the coefficients will be printed 

the means and standard deviations unsealed. A table showing 

the number of positive responses for each item will be printed. 

To read data in seconds; number of persons (.007 var. + 020 ) 

Calculate and punch in seconds; (.055 + «017 d n + .0005 number 

of persons) per coefficient. 
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Add 2 minutes for miscellaneous reading and punching. 
FF STOPS FF location ONI - character read not .0, 1, N, F or L. Item 

and person numbers printed* Raise white switch up and resume 
checking of data. 

FF location ONJ - a person's number of items does not agree 
with row 1. Computer will print item and person number. 
FF location 1F8 - sum check failure on master readin. 
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KSL 2,50 
FISHER'S Z TRANSFORMATION t (Z - l/2 en(l + m)/(l - m)) 

PURPOSE This program allows you to determine Fisher's Z which can 

be used for averaging coefficients of correlating and in 
testing the significance of differences between correlations . 

OPERATOR'S INSTRUCTIONS Stops 

1, Master KSL 2.^0 2kOQ6 

2. Data tape (black switch up for each 

row) 2k006 

As many data tapes can be read in without rereading master. 
When machine reads ah N-— a 2k006 stop occurs. When a J 
is encountered, an OF occurs. On this OF more data tapes 
can be read by raising white switch up and down. 
PREPARATION OF TAPES 1. Master - duplicate library copy KSL 2.50 

2. Data tape - the data tape consists of a series of signed 
fractions scaled 10 and terminated by an N or NJ (see 
Operator ' s Instructions ) • The typical data tape is the 
output of K-8 or K-9 or any other correlational routine. 
The data tape can be triangular output (insert an N) 
on rectangular output (N's already present). 
TIME ESTIMATE About .415 seconds for each number read in. 
FF STOPS FF location 1L5 - sum check failure on read in of master. 



KSL 2.51 
COMPUTES A TRIANGULAR CORRELATION MATRIX CONTAINING VARIANCES-COVARIANCES 

PURPOSE Given a triangular variance-eovariahce matrix this routine 

computes and outputs a triangular correlation matrix where 

each correlation obtained is = — «- where (T . . is covanance 

^i^j 1J 

of i and jo <j- . is standard deviation of i and (T . is 
standard deviation of j . The number of columns must be < 40 . 
OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 2.51 2 ^ 001 

2, Parameter tape 2403K 
5. Data tape (black switch down) OF 

PREPARATION OF TAPES 1* Master - duplicate library copy KSL 2.51 (Be 44) 

2 . Parameter - + c + d n N where c = number of columns must 
be < 40 d is the number of decimal places in output 
must be < 12. 
3„ Data tape .-■ a. triangular variance-covariance matrix of 

the form a~y,, <r ~\2 9 <r ~22' <I ~23' * ° " ' ^nn^ ' 
Samples can follow one another. The last sample should be 

followed by an NJ. 

No time estimate or FF stops given . 



KSL 2.6l 
CHI-SQUARE FOR M BY C FREQUENCY TABLES 

PURPOSE This routine will calculate a set of chi-square values for 

a set of cell frequencies for a table of M rows and C columns 
If desired, expected frequencies can be obtained . These are 
found by multiplying the ith row sum by the jth column sum 
and dividing this product by the total frequency of the table 
For each table a chi-square value, the degrees of freedom 
will be punched on tape, as well as the expected frequencies . 
(See Operator's Instructions). 2rc + r + c < 680 is the 
capacity. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL 2.6l 2k08K 

2. Data tape OF 

To get expected values, chi-squares and degrees of freedom 
white switch up and down to read in data. To get only chi- 
square degrees of freedom read data in with black switch. 
Additional data tapes can be read in by raising white switch 
up and down and then using either white switch or black 
switch as indicated above. 

PREPARATION OF TAPES 1. Master - duplicate library copy KSL 2.6l 

2. Data tape - the cell frequencies are punched on tape 
as positive signed integers. Each row is terminated 
by an N, and the last N of the table is terminated by 
a J. After the last table on the data tape, an L is 
placed. 

FOR EXAMPLE +3 + .+29 +4N 

+4 +52 +1N 
+16 +2 +9N 
+111 +49 +1NJ ..... NJL 

TIME ESTIMATE Chi-square and degrees of freedom in seconds per table 

(.033r c + 1.2) 

With expected values,,: -chi-square and degrees of freedom, 
in seconds per table (,210r c + 1.8). 



FF STOPS FF location 125 - sum check failure on master readin. 

2009K - number of elements in row does not agree with number 

of elements in first row. Machine prints row number, 

2*K)8S - If frequency is negative due to data tape error - 

machine prints row number. 

2koQS - If chi-square value is undeterminate (row or column 

sum totals zero) machine prints error undeterminate. 

To continue reading remainder of data with the three listed 

black stop error j raise black switch. 
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PURPOSE 



KSL 2.70 
D. STATISTIC 

The program computes the index 






n (a . - a . Y 

1 x vi vj ' 



1/2 



OPERATOR'S INSTRUCTIONS 

1 
2 
3 



PREPARATION OF TAPES 



EXAMPLE 



for all i and j with j < 1 = 2 3 . .., m. The results are 

printed in a triangular form where columns are associated 

with index j and rows with index i. This set of indices 

is computed successively for each of a sequence of persons. 

Size of each item must be a single digit. Number of items 

per description < 100 and number of descriptions < 12. 

Stops 

Master tape KSL 2.70 2k0k6 

Parameter tape 3404L 

Data (black switch down) read in order 24o4K 

if black switch up 

Master ■- duplicate library copy KSL 2„70 

Parameter tape - the number of descriptions per person 

(m) and the number of items per description (n) are 



1. 
2. 



given in this form 



+00M 
+00N 



For example,, 10 



m; 99 = n would be +010 

+099 



3 o Data tape - all the item scores for one description 
given by a person followed by all the item scores on 
next description by that person, and so on until all the 
data for that person is given. An N is punched at the 
end of each person. An NJ is punched after the last 
person on the tape. 

Consider a problem in which there are two persons 3 each 

filling in h items on each of five descriptions . Then the 

parameter tape would bes 

space +005+004 spaces and the data tape would bes 
d l d 2 d 3 d 4 d 5 

P, 1234 1234 6^4l 3256 4711 N 



P 6666 1111 2244 4421 1244 NJ 



The results would "be printed as follows 
001 = person number 

m m+1 m+2 m+3 
m+1 00001 
m+2 06633 06633 
m+3 03464 03464 06633 
m+4 O6856 06856 04123 O8185 

002 

10000 

06325 04472 

07000 04359 0^583 

07000 04359 01000 05099 

„ ^ * m(m-l) N 
TIME* ESTIMATE Number of seconds = — ^ x — 

where m = number of descriptions and 
n = number of persons 
FF STOPS None listed 



PURPOSE 



KSL 2.9^ 
AGREEMENT PATTERNS,, CODE H 

This program finds the highest agreement score among all 
the persons. An agreement score is a tally of the number 
of items on which two persons agree. Additional persons 
are selected as long as the joint scores multiplied by the 
number of persons in the group is increased * When the 
score no longer increases., the persons and joint scores 
(agreement score) are printed, and the joint pattern of 
responses is removed from the measurement matrix forming 
a residual measurement pattern matrix o The process con- 
tinues by again finding the highest agreement score in 
the residual matrix, repeating the above steps resulting 
in a second residual matrix. Whenever the agreement scores 
become lower than some predetermined value, the process 
is at an end. The capacity can be determined by the 
following; 

780 > ~r P + P (P > 25)| 758 > Yo p + p ( p < 2 5) 



19 _ . _ % . - -„, ,.- 19 

where I = items and 
P = persons. 
OPERATOR'S INSTRUCTIONS 

1. Master tape KSL 2,$k (part i) 

2. Data tape 

3« Master tape (part II ) resume where 



Stops 
2^013 
243F7 



stopped on (l) above 2^5*7 

Black switch up to begin calculating. If the amount of 
time required to reduce the matrix to zeros is going to 
be prohibitive, insert the following parameter with a 
black switch up and proceed by black switching to begin 
calculating s 

0084+ 

00F OOSF 

2685N where S is the agreement 
score to stop at. 



PREPARATION OF TAPES 



FORM OF OUTPUT 



TIME ESTIMATE 



1. Master - duplicate library copy of KSL 2«9^ 

2. Data tape - item responses can be 1, 2 or 0. No signs, 
a - at end of each person and an N, J, F or L at end 

of the measurement matrix, 
3o Parameter tape - see OPERATOR'S INSTRUCTIONS. 
An item count is made on each person. If data tape is 
correct (only 0, 1, or 2's and the same number of items for 
each person) , the machine will print number of items and 
persons before beginning punching of agreement scores. 
If there are errors, the machine "will punch the type of 
error and the person in which the error occurred. The 

machine will then stop on a 24014. 

PI 
To read in data tape: ^ro seConds 



FF STOPS 



Calculation: 

Printout; 
None 



2P + 221 seconds (i = number of items 

P = number of persons) 



- PI 

= •=■?? seconds 
50 



KSL 4.00 
FREQUENCY DISTRIBUTIONS 



PURPOSE The routine will read a set of n signed integers . These 

integers will be grouped in ranked categories where the 
lowest boundary of all categories is set by the user . 
When the letter N is read at the end of a signed integer , 
the routine will print a frequency chart with a bar graph 
alongside. The width of the bar graph is automatically 
adjusted so as not to exceed the width of the teletype 
paper. The maximum number of categories is 100. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL ^.00 3^0L7 

2 o Parameter tape 2411N 

3. Data tape 2^-111 

The machine will stop on a 2^11N at the end of each group 
to be categorized o A black switch up will read a new group 
using the same parameter. If a series of n integers are 
to be categorized using the same lower boundary and in- 
crement between categories, a black switch can be used. 
However, if a new parameter is to be used, white switch 
up and down to read in the new parameter at a 2411N stop 
or when the machine reads a J, at which time it will stop 
on an OF. 

1. Master"- duplicate library copy KSL ^.00 

2. Parameter tape - the tape consists of two signed integers 
and is terminated by an U. + B + I No J: ^ is ^ e 
lower boundary of the lowest category. I is the increment 
between boundaries . 

100 integers ranging -9 and +90 in value. Parameter might 
be -9 + IN. This will give a frequency distribution over 
99 categories. 1000 integers ranging between -100 and +200 s 
Parameter might be -100 + ION. This would give you frequency 
distributions over the intervals -100; -90| -80| -70| •••;> 
+180| +190 | +200 - a total of 31 intervals with actual 



TAPE PREPARATION 



EXAMPLE 



TIME ESTIMATE 



numbers grouped to fall within these respective categories . 

If you should fail to list the lowest number in the n 

integers j the numbers falling below the lowest category 

given in the parameter and above the range of 991 will be 

printed out as an ERROR and the remaining numbers will be 

categorized. 

Data tape - the data tape consists of signed integers 

terminated by an R. 

[(4d. + ho) n + 650 c] 10"^ seconds 

where 



d. = digits per integer 
n = number of integers 
c = number of categories 
FF STOPS FF location 1V3 - sum check failure on master read in .- read 

again. 
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KSL ^ o 05 
RANKED NUMBERS 

PURPOSE This routine takes in a series of numbers and ranks them 

according to a preset parameter. Each row is ranked and on 
the output is designated as a set. Two columns appear on 
the output tape-o One contains the row numbers of the variable 
and the positive numbers ranked in descending order. The 
left column contains the negative numbers in descending order. 
The problem is limited to 26 columns of 200 numbers each. 
Output is in four digit form. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 4.05 243F7 

2 . Parameter 20019 

3. Data OF if J at end of tape 
A new parameter and data tape may be read in at this point 

by rereading master. 
PREPARATION OF TAPES 1. Master - duplicate library copy KSL 4.05 

2. Parameter - 0012+ 

OOF OOS^F 

OOF OOSgF 

OOF OONF 

OOF OOP^F 

2620N 

+ L A N 

where S^ = begin set (usually l) 
-p 
S = set end. (S fi and S is affect give the limits 

of number of columns) i.e., S = number of columns 

N = number of signed elements in each column 

P = limit on number of numbers wished printed 



(P < N). [If you have 200 elements (N) and 
x n — 

wish to know highest 30 elements (both + and ^) 



set P =30] 

n ^ 



L * = Lover limit (absolute value) cut off point. This is 
obeyed before P is obeyed. For example, you have a 
column of 200 numbers. You are interested in numbers 
with an absolute value of 0200 and over. When all 
numbers of 0200 and over are printed, regardless what 

P is. machine will go on to next column. However. 

n ' ■ 3 

if you ask for a P N of 10 and there are 20 numbers 
of absolute value over 0200, it will only print out 
the highest 10. 
*N0TE (if your data is scaled 10 , 10 , etc • , take this into 

account when setting L. - i.e., a L. of 2000 with a column 
of numbers scaled 10 - K-8 output .- will cause no numbers 
to be ranked - they all lay below this cutoff point). 
3» Data - this is a series of signed numbers, each row or 
column terminated by an N. They may be scaled or not.* 
TIME ESTIMATE No time estimate given. Generally this is read in time 150 

characters per second and print out time 50 characteis _ per 
second. At one/half minute for each set for calculation. A 
six digit, 100 number, 20 column problem would be approximate 
(printing out only top 20 numbers). 

V s Nx S E , T , 6 x 100 x 20 o_ 

— _ - . ^ — - seconds = Input =* t~ — _ ■ ■ . - - ■ = oO seconds 

150 * 150 

— r— minutes = Calculation = — -j- = 5 minutes 

T- X P X Q_ I, _ _ _ _ 

n E , _ . , k x 20 x 20 , lQ 
— " n ■■ seconds = Output*-* — — ~-^- = 40 seconds 

Total time = 7 minutes 
FF STOPS No FF stops given. 
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KSL 4,10 
STANDARD SCORES 

PURPOSE This routine takes in a raw score matrix by row and 

standardizes these scores by columns according to preset 

parameters provided., The standard score obtained is 

— ' 

^pr- (x ~ m) + m ! where G~ '■ is standard deviation of the 

standard s cores S(T - is standard deviation of the scores, 

x is the raw score, m the mean of raw scores and M' the 

mean of standard scores. M' and<3~" ' are preset by user. 

The number of v (elements in a row) < 1^5. Sample size is 

unlimited, the output is signed, *** each row terminated by 

an N and the last N followed by a J. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master tape KSL 4.10 24001 

2. Parameter tape 24053 
3« m -CTtape 3406k 
4. Raw score tape 24o6^ 

(a) one more person's raw scores may be read in by raising 
black switch. 

(b) to read in a new parameter, m -CF" tape and raw score 
tape, place parameter tape in reader, raise white switch up 
and down, then repeat steps 3 and 4. 

PREPARATION OF TAPES Master - duplicate library copy KSL 4.10 

Parameter - d n space n space +m' +** <T"* ' ** H .or J 

where d n is the number of decimal places in output < 11. n 

is the number of rows, people or sample size, m' is the mean 

of the standard scores . g— - ' is the standard deviation of 

the standard scores • 

If you wish the output in columnar form use an N; if you 

wish it printed in rows use a J. (There are no carriage 

returns provided in output by rows except at the end of a rov$. 

M -*—tape - this tape contains the means and standard 

deviations of the raw scores. These are output of K-8 or 

KSL 2.03 with an N added after the last cr~o 



If these are punched directly on tape by the user, they 
must be signed and follow this pattern 

M_ o~ of 1st variance 

yL q~ of 2nd variance 

M-, <r~- of 3rd variance 
3 3 



ESTIMATED TIME 

FF STOPS 
NOTES 



M <3~ of nth variance 
v _y 

An N must terminated the M - 0~" tape . 

The number of digits in each number is < 11. 

Raw score tape - this tape is punched person by person . (As 

for K-8 data tape). Each element (variable) must be signed 

and have < 12 digits. An N must terminate each row.* The 

same number of variables must be in each row Scaling is 

unimportant by row. It is important by variable. For 

example s there follows a sample 3 variable k person matrix s 

+001 =02 +Q011N Not scaled by rows 

+010 -=21 -0101N Scaled by columns 

+111 +00 +1246N 

+099 -32 



Approximately 



0321N 
16 



100 



seconds per raw score 



FF001 means wrong number of variables in a row. 
* If an F terminates a person scores (in place of an N), 
the computer will stop^ allowing another data tape to be 
read in. Raise black switch to begin reading in next tape. 
** Suppose M' is set equal to +50.00 and o" * = +10.00. 
These will be punched on the parameter tape as +5 and +1. 
Likewise 9 a M' = +.5 and o~~ ' = »1 are punched as +5 and 
+1. However 3 if M 1 = +5.00 and q- ! +10. 00 , these are 
punched as +15 and +10. 

*** If the raw scores are scaled by 10~ x ,, then the standard 

~x+l 
scores will be scaled by 10 



PURPOSE 



KSL 4.15 
CUMULATIVE PROPORTION 

This routine reads in a matrix of raw scores, and converts 
each element to a cumulative proportion value lying between 
0*000 and 1.000. The user has two options on output . 
(a) The cumulative proportions may he printed out in the 
form of a frequency table giving score, frequency and cumulative 
proportion table or (b) the score matrix may be punched out, 
each element of the raw score matrix being replaced by its 
corresponding cumulative proportion value. If B option is 
used, the results from this program may be fed directly into 
KSL 4.l6 or KSL 4.17° If option A is used, no further use 
can be made on the computer with the output tape. The 
number of scores that can be ranked at one time is 650 or 
lesso The number of groups that may be run is limited by 
only the amount of time allowed on Illiac for one run. The 
output consists of signed (+) four digit numbers, each row 
terminated by an N, the final row terminated by an NJ. 

OPERATOR'S INSTRUCTIONS §i°£S 

1. Master KSL 4.15 2 ^°5 J 

2. Parameter 2*1-013 

3. Data 243F7 

A new parameter may be read in with a black switch up . 
Then, a new data may be read in. 
TAPE PREPARATION 1« Master - duplicate library copy KSL 4.15 

2o Parameter - 006-4- 

00F 00 cF 

OOF 00(*)F 

OOF OOrF 

2419N 

where 

c = number of elements in one row or columns in raw score 

matrix. 
* = ty p e f output 1 a will give cumulative scores replacing 
raw scores. This output usable in KSL 4.l6 or 4,17$ or 
1 will give cumulative scores in a frequency table. 



KSL 4.15 
CUMULATIVE PROPORTION 

PURPOSE This routine reads in a matrix of raw scores, and converts 

each element to a cumulative proportion value ;Lying between 
0.000 and 1.000. The user has two options on output . 
(a) The cumulative proportions may be printed out in the 
form of a frequency table giving score, frequency and cumulative 
proportion table or (b) the score matrix may be punched out, 
each element of the raw score matrix being replaced by its 
corresponding cumulative proportion value* If B option is 
used, the results from this program may be fed directly into 
KSL 4.l6 or KSL 4.17. If option A is used, no further use 
can be made on the computer with the output tape. The 
number of scores that can be ranked at one time is 650 or 
lesso The number of groups that may be run is limited by 
only the amount of time allowed on Illiac for one run. The 
output consists of signed (+) four digit numbers, each row 
terminated by an N, the final row terminated by an M. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 4.15 2405J 

2. Parameter 24013 

3. Data 243F7 

A new parameter may be read in with a black switch up. 
Then, a new data may be read in. 
TAPE PREPARATION 1. Master - duplicate library copy KSL 4.15 

2o Parameter - 006+ 

OOF OOcF 
OOF 00(*)F 
OOF OOrF 
2419N 
where 
c = number of elements in one row or columns in raw score 

matrix. 
* = "type of output^ a will give cumulative scores replacing 
raw scores. This output usable in KSL 4.l6 or 4pl7$ or 
1 will give cumulative scores in a frequency table. 



3o Data tape - this tape consists of a raw score matrix „ 
The usual raw score matrix has £ columns and r rows where 
the £ are variables and the r are persons o IF YOUR RAW 
SCORE MATRIX IS IN THIS FORM YOU MUST TRANSPOSE IT BY KSL 5*30, 
This will give you a raw score matrix of £* columns and r* 
rows where the c* are now people and the r* are variables » 
Therefore,, the cumulative scores will be given over all 
persons for each variable <> Each element must be preceded 
by a plus sign., Each group of scores to be ranked must be 
terminated by an No If your raw score matrix has values that 
are negative 3 add a constant to each element so that all 
elements are positive <> 

TIME ESTIMATE (5 x r x c x highest score x 10~^) seconds to calculate. 

{See note) 
(ol x r x c) seconds 

FF STOPS FF015 - number of elements in row does not agree with number 

in first row e 
FF06L - a negative number - white switch up and down - 

minus sign will be replaced by a plus and problem 
will continue o 

NOTE Since this program gives cumulative proportions according 

to the position of the score in the group, a very high 
value may be replaced by a smaller one - thus saving much 
machine timeo For example,, scores are +1+3 +0 +2 +999 • 
The 999 will receive the highest cumulative proportion value 
as long as it is higher in value then the other k- scores „ 
Therefore 9 a 9 maybe used to replace this value without 
affecting the final results , but with the advantage of 
saving machine time 



KSL 4.l6 
T-SGORE PROGRAM 



PURPOSE This program converts a set of signed (+) four-digit 

cumulative proportion (results of KSL ^.15) numbers to 
normalized T-scores. It is essentially a table-look 
up operation . The output tape consists of signed (+) three 
digit numbers, ranging between the values of 19.0 and 80.0. 
Groups of scores (n people over 1 variable) are terminated 
by an N. The last row in the matrix is terminated by an 
NJ. This output may be transposed by KSL 5. 30 and is then 
ready for use in K-8. No lltmit to size of matrix. 

OPERATOR'S INSTRUCTIONS Stops 

1. Master KSL 4„l6 2^015 

2 . Data OF 

1. Master - duplicate library copy KSL 4.l6 

2. Data - normally this will be output from KSL 4.15. 
However, it may be prepared on the teletype equipment. 
If so, it consists of rows of signed (+) four-digit numbers 
lying between +0.000 and +1.000. Each row is terminated 
by an N, the last row is terminated by an NJ. 
Approximately 1/3 second per T-score (r x c). 
FF stop - a + sign is missing. 

An N on data tape causes N to be printed on output tape. 
A J on data tape causes machine to print a J on output 
tape, then come to an OF stop. 



TAPE PREPARATION 



TIME ESTIMATE 
FF STOPS 



KSL 4.17 
STANDARDIZED T-SCORES 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



TIME ESTIMATE 
FF STOPS 



This routine takes signed (+) four-digit cumulative 
proportion numbers (results of KSL 4.15) and converts 
them to standardized scores - punching them out on 
tape. The standardized scores lie between the values 
of -3.00 and +3.00. The output tape consists of signed 
three digit numbers . Each row is terminated by an N»- 
The final row is terminated by an NJ. There is no limit 
to the size of the matrix. used. The output may be trans- 
posed by KSL 5.30 and is then ready for K-8. 

Stops 

1. Master KSL 4.17 2*1-015 

2. Data OF 

1. Master - duplicate library copy KSL 4.17 

2. Data - normally this will be output from KSL 4..15. 
However _, it may be prepared on the teletype equipment. 
If so, it consists of rows of signed (+) four-digit 
numbers lying between +0.000 and +1.000. Each row is 
terminated by an N, the last row is terminated by an NJ. 
Approximately 1/3 second per T- score (r x c). 

FF stop - a + sign is missing. 

An N on data tape causes N to be printed on output tape. 
A J on data tape causes machine to print a J on output 
tape, then come to an OF stop. 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



PARAMETER TAPE 



KSL br.50 
LIMITED INFORMATION ESTIMATION, SINGLE EQUATION 

This routine estimates parameters in economic models "by 
the limited information single equation method » Consider 
the equations y = B^ y 2 + B^, y, + ... + c^ z^ . .. + u^. 
The routine estimates the B's and C's and also calculates 
their standard error, A parameter tape specifies which 
of the endogenous and exogenous variables are to be in- 
cluded in the equation under study. To find estimates for 
other structural equations on the same data, use a different 
parameter tape.. 

Stops 
lo Master tape KSL 4.50 3^3KS 

2. -Data tape I 2*1-108 

3a. Parameter 2^108 

New parameters for Data I can be read in with black switch . 
To read in a new data tape II, move white switch up and 
down then black switch in new parameters. 

1. .Master - duplicate library copy KSL 4.50 

2. Data - the data tape consists of one or more covariance 
matrices . (The order or -size of these covariance 
matrices must be < 14). Each covariance matrix in the 
set must have the same number of elements. The elements 
are punched as signed fractions, and the matrix is 
punched in triangular form (such as the output from K-8) 

The final matrix in the set S is terminated by an NJ. The 

largest number of matrices in any set S, must satisfy this 

o 
in equality s S(n + n) < 18,600 where n is the number 

of rows in one matrix. 

Each parameter tape determines a structural equation and 

operates in turn on each of the covariance matrices in the 

set. The parameters for different equations in the same 

mode may follow one another on the same tape. 



+a +a p + ...+aJI where a f s represent the row numbers of 
the covariance matrices which are to "be 
the endogenous variable in the equation. 

+b +b +...+b N where b's represent the row numbers of 
1 2 g 

the covariance matrices which are to be 

the endogenous variables in the system. 
+c +c p +...+c h N where c's represent the row numbers of 
the covariance matrices Which are to be 
the exogenous variables in the equations. 
. +tN t is the number of observations used. 

Limits of f, -g, h, and t are 
2 <£ <6 1 + -f + g < t 

2 < g < 8 f + g < N. (size of matrix) 
0<h<6 
TIME ESTIMATE Read master - 90 seconds 

Read set of covariance matrix = 5 seconds/matrix 
Calculate and punch results = 15 seconds/matrix. 
FF STOPS FF location 3KK .- sum check failure on master read in. 

Symbols may be punched in the results . The machine does 
not stop but reads in next matrix. The meaning of these 
are as follows? 

0/xxx number of elements in a covariance matrix disagree 
with number in first matrix of the set. The matrix 
number 9 xxx^ is in sexadecimal, 
l/ a submatrix is singular or nearly so. 
,2/ overflow on an element in matrix ..multiplication. 
.3/ failure in rescaling VL R or R~ W. Scale the 

covariance matrix down and try again. 
k/ failure in adjusting estimates or standard errors 

to proper scaling for printing. .Try a correlation 
matrix instead of a covariance matrix (see complete 
library writeup for description on how to do this). 
.5/ failure in forming F/x" 1 . The value of \ .is pro- 
bably near zero. 



- 2 - 



6/ trial vector of (iT ¥ - XI ) u = has not 

converged after 30 iterations » 
7/ failure in rescaling elements of F, %. 
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PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



KSL 5.02 
MULTIPLICATION OF A MATRIX BY ITS TRANSPOSE 



T T 
Given a matrix A. ., this routine will form AA or A A 

T ^ 

where A is the transpose of the matrix. If the columns 

in the original matrix A. . are < 6k. then the number of 

rows must "be < l60; if the columns are < 128 , the number 

T 
of rows must be < 80, The output matrix of either AA 

m 

or A A will be symmetric. (Lower triangular half = upper 
triangular half). 

Stops 

1. Master tape KSL 5-02 3^087 

2. Parameter tape 2^0K7 

3. Data (Matrix A) tape 2^087 

At stop 2^087 , a new parameter tape can be read by raising 

the black switch. Follow this by reading in a new data 

tape. 

1, Master - duplicate library copy KSL 5 .02 



2, Parameter - d^ JtF x L + sS 



V 





number of decimal digits desired in output 

T 
t : t = 0, A A will be formed 

T 
t = 1, AA will be .-formed 

x : x determines the form of output matrix 

x = output in triangular form by columns with 

carriage returns after each element, 
x = 1 output in triangular form by rows punched 

across the page* 
x = 2 output will be square by columns with a 

carriage return after each element and an 

N after each column, 
x = 3 output will be square by rows across the 

page with an N after each row. 
+s: +s or -s is the integer exponent of ten. To 

scale by 100 +s = +2; to scale by .01, -s = -2. 
Remember, matrix elements read in to machine are 
treated as if decimal point is after sign. (See 
FF stops ) . 



TIME ESTIMATE 



FF STOPS 



+Oss If no scaling is desired , sign precedes the 0. 
3. Data tape - the elements of matrix A are punched by 
rows as signed fractions with an N at the end of each 
row. The final N must be followed by a J. If an 
F is punched instead of an N, the computer will stop 
on either 300K8 or 300SS. By raising black switch 
another tape containing the remainder of matrix A 

may be reado 

T 
AA square or triangular form in seconds? 

p 

(number of rows) (.02d + .002 number of columns + .12) 

+ number of rows (»005d n number of columns + .018 number 

of columns + .03) +^0. 

T 
A A in square or triangular form in seconds: 

substitute number of columns for number of rows and number 

of rows for number of columns in formula above. 

FF001 location OSF - number of elements in rows are not equal 



FF002 location 11F, 171 
FF003 location 112, 175 



FF004 location 13K 
FF005 location 13J 

FF006 location 3F5 



overflow in an element of matrix c 
To read in a new scaling parameter 
-sS only, move white switch up and 
down j the white switch again to 
eliminate having to reread data* 

Output will ensue » 

T 
sum check failure in matrix A or A . 

white switch three times 5 re-insert 

data, and then black switch to start. 

sum check failure on master read in. 



KSL .5,11 
MATRIX INVERSION 

PURPOSE This routine takes either a triangular matrix (output 

from K2, K8, K9) or a square matrix and computes and 
punches out the diagonals of the inverse matrix and/or the 
inverse matrix o The inverse matrix can "be punched out 
in the form of a triangular and/ or square matrix « 

OPERATOR'S INSTRUCTIONS Stops 

1. .Master KSL 5-H 3^09S 

2. Parameter 200K3 

3. Data 2^09S 

To run new data tape on same parameter, white switch up 
and down. 
TAPE PREPARATION 1. Master - duplicate library copy KSL 5.11 

2, Parameter d~ space N space X space 

where d n = number of digits in each element of output 
N = order (number of rows or columns of matrix) 
X = - diagonals only punched 

1 - diagonals and triangular output 

2 - diagonals and square output 

3 - diagonals and triangular and square output 
(The above are used if a triangular matrix is data matrix). 

X = 10 - diagonals only punched 

11 - diagonals and triangular output 

12 - diagonals and square output 

13 - diagonals and triangular and square output 
(The above are used if a square matrix is data matrix) « 

3. Data tape - the data tape can be either a triangular 
or square matrix « If it is triangular the termination 
symbol at end of entire matrix must be an N, J, F, L 

+ or -• If it is a square matrix, each row must be 
terminated by an N, J } F, or L. In either case, each 
element must have a sign and not more than 12 digits 
in an element , 



TIME ESTIMATE Not given * 

FF STOPS FF018 - diagonal entry negative "j vmic h error 

FF01N - diagonal is zero J 

FF01+ - singular matrix 

FF019 - any diagonal is negative - to "be used as pivit„ 
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PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



KSL 5.16 
RESCALE ML3 OUTPUT 

This routine takes the output of ML3* scales all elements 
by the same power of 10 and punches the data, column by 
column, as a square matrix with N's after each column and 
a J at the end of matrix. 

Stops 

1. Master tape KSL 5.16 2k3Ff 

2. Parameter tape 2^06f 

3. Data tape 2^3F7 

A new parameter and data tape may be read in with a 
black switch. 

1. Master tape - duplicate library copy KSL 5»l6 

2. Parameter - 00*H- 

OOF 00 (d . )F 
OOF 00(d Q )F 
OOF 00(c)F 
OOF 00(r)F 
OOF 00(p)F 
26100N 

+ (s) 

where d, = number of digits in each element input 
d~ = number of digits desired in output 
c = number of columns of each matrix on problem 
tape (if same parameter is used for more 
than one problem - the number of columns and 
rows must be identical for each matrix), 
r = number of rows of matrix 
p = number of matrices on problem tape 
s = the scaling of Ml 3 output (given a 10 column 
matrix from ML3, the scalers at the end of 
each column may vary between +1000 and +00000 
... 1. If for the matrix the smallest scaler 
is +00001, then (s) would be this number). 
3» Data tape - this is simply the output from ML3» 



TIME ESTIMATE 



FF STOPS 

TOTtfmv 



Essentially input-output time. 

c x r x d, 

2 

Input time in seconds = t-Tq — — 

c x r x d 



Output time in seconds = 



"0 



50 



where c, r } d, and d are defined under parameter's above 

None except in library subroutines • 

The scaling factor of this output can be determined by 

the formula io^ 7 " 1 where K is determined by the scaling 

of input into Ml 3 - (l0~ ) and y is the scale factor of 

output of ML3» 

■ ■ : -1 
For example: input of Ml 3 scaled; 10 

scale factor (+00l); 10" 5 

KSL 5«l6 output scaled 10 = 10 • 
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KSL 5.20 
MATRIX ADDITION OR SUBTRACTION 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



Given two matrices A. . and B. .. equal in number of rows 

and columns , this routine will add the elements of A . . 

to B. . and print out the new matrix C„ . . The routine 

U U 

will also subtract the elements of B. . from the elements 

ij 
of A. . and print out the new matrix D„ . . It is possible 

to sum a series of B. . matrices before subtracting or 

ij 

adding the elements of B. . from or to A... The A. . and 

U ij ij 

B. . matrices can be rescaled by any power of 10 up to 

10 to prevent overflow of the elements of C. The size 

of the matrices must satisfy this inequality i(j-f-l) < 10^2^0 

where i = number of rows and j = number of columns in one 

matrix which must be less than 400. 

Stops 

1. Master KSL 5°20 3^022 

2. Parameter tape 2403N 

3. Data tape A 3^057 



ka. Data tape B. 

hh • Data tape B,. 

Data tape B. 



1 



n 



24057 
24057 
0F000 



To begin a new problem at stop 0F_, move the white switch 
up and down. The computer stops on 2^022 . Insert the 
new parameter tape and raise the black switch. 

1. Master - duplicate library copy KSL 5«20 

2. Parameter tape - this tape consists of three signed 
integers terminated by an No They are in the following 



order : 



■KL + n + sN where 



d n = number of decimal places in output 
n = number of B matrices . A positive sign indicates 
additiono A negative sign indicates subtraction. 

s = the scaling integer. The A and B matrices are 

+ 

multiplied by 10— s during input. If no rescaling is 

desired (s=0) 5 this parameter may be omitted. If 



both n and s are omitted 3 two matrices 5 A. . and B 



U 



ij 



will be added without rescaling, 



3o Data tape - the elements of A. . and B. . matrices are 
punched by rows as signed fractions with 12 or less 
digits in each element <> Each row is terminated by an 
N and each matrix is terminated by a Jo If no F is 
punched in place of an T$ 9 the computer stops j, allowing 
a new tape to be read in. 
TIME ESTIMATE ,005n(d + 2)ij + .017 (d + l)ij seconds . 

nii.ere n = numoer Ox a ma ^nces 

d = number of digits per element in decimal output 
i = number of rows in matrix and 
j = number of columns in matrix 
FF STOPS FF000 location 0k6 - elements of A vector overflows into B 

vector | white switch to read another 
row of A. 
FF001 location O56 - number of elements of rows of A not 

equal . 
FF002 location 065 - number of elements of columns of A 

and B not equal o 
FF003 location 071 - number of elements of rows of A 

and B not equal. 
FFOO^- location 137 - master tape sum cheek failure 3 reread 

master tape. 



KSL 5<30 - 259 
MATRIX TRANSPOSITION AND RESCALING 



PURPOSE 



OPERATORS INSTRUCTIONS 



PREPARATION OF TAPES 



This routine reads in a rectangular matrix and punches 

it out transposed and/or rescaled according to user's 

-2 -1 
designation. The elements may "be rescaled 10 ,10 , 

10 +1 or 10 . The form of output on teletype page print- 
out may be columnar, row or in a row with a space separating 
elements. It does not provide carriage returns at the 
end of a row of teletype print except after a complete row. 
That is - if the total number of characters (sign included) 
on a line or row of teletype print exceeds 70, the resulting 
2nd, 3rd etc., lines left-hand margins will be ragged in 
printing. Neither R (rows) nor C (columns) can exceed 775 
elements* R x C < 102^4-0, This routine is especially 
helpful in going from the output of one routine to input 
for another. (For example, KSL 1.20 to KSL I.96). 

Stops 

1, Master tape KSL 5-30 3^0FN 

2, Parameter tape 2^0F8 

3, Data tape 2*K)F8 

An additional data tape may be read in with a black switch 
at this point. It uses the same parameter. To change 
parameters at stop on step 3, insert parameter, raise white 
switch up and down and repeat step 3« 

1. Master - duplicate library copy KSL 5*30 

2 . Parameter 



consists of d n space X S P 



. „ decimal places wished in output, (if you are reading 

+2 
in elements with 6 decimal places and scaling 10 you 

-2 
would ask for \ decimal places; if you are scaling 10 ask 

for 8 decimal places; if you are just transposing, ask for 

6 decimal places ) , 

x =i to transpose matrix; x = 1 to rescale only and 
x = 2 if the matrix is to be transposed and rescaled. 



TIME ESTIMATE 



FF STOPS 



NOTE 



EXAMPLE 



If x = 0,, omit S 

If s = +2 (sign Important) elements will be scaled 10 

.+1 



+2 



s = +1 elements will be scaled 10 
s = -1 elements will be sealed 10 



-1 
,-2 



PS 



s - -2 elements will be scaled 10 
If P = E 3 output appears printed by columns 



P = 



output appears printed by rows 



P = F, output appears printed by rows with a space 
separating elements . 
Data tape - the data tape must have an N terminating each 
column or row A J follows the final N. Each element is 
a signed fraction < 12 digits » For transposing, the number 
of elements in each row must be equal.. For scaling alone, 
this is not important . A large tape may be broken Into 
smaller ones at the end of a row or column and an F sub- 
stituted for the No The computer stops and a second tape 
may be read in with a black switch. 

In seconds = .25 R x C where R = number of rows and C is 
the number of columns o 

4F location 053 - number of elements in row ^o not agree 
with number in first row. 

FF location 100 - sum cheek failure master read in* 
If you wish to scale a triangular matrix (K-8 output), 
place an NT at end of matrix and use KSL 5.50. Do not 
try to transpose anything but a rectangular matrix. 
52 -U Parameter - transpose and scale 10™ j print by rows 

to 5 places, [i.e., one leading zero plus 4 digits 

in column 10] . 

DATA INPUT 

+3-4+06+100-29+30+1+2-2 39+0632N 

+I+5+IO+099-3O+31+2-I-225+HIIN 

+2+2+ll+098=32+29+l+2+321=2103N 

-3+0+07-001-45+52+2-2-112+12 46N 

-2-2-2I+056+01+31+3+3+099-0029N 

+4-1+42+113-18+32+1+4-021+0632N 

+9-^-2l+2l0-l9+25+i+2+l25+2567N 
-3-9-32+098-32+42+1-2+625-52^ 
+9-9+61+217-41+32 -1+6 -2 35+1112N 

+8+2+07-304+29f55=l+l+325+H13N 
+7+4+09-002+42+24-2+2-249-2468N 

T 



OUTPUT RESULTS 

+03000+01000+02000-03000-02000+04000+0 9000-03000+09000+08000+07000N 
-04000+03000+02000+00000-02000-01000-04000-09000»09000+02000+04000N 
+OO6OO+OIOOO+OIIOO+OO7OO-O2IOO+042OO-O2IOO-.O52OO+O6IOO+OO7OO+OO9OON 
+01000+00990+00980-00010+00560+01130+02100+00980+02 71O-03O4O-OOO2ON 
-029OO-O3OOO-O32OO-O43OO+OOlOO-Ol8OO-O19OO-O32OO-04lOO+O2 900+OteOON 
+O3OOO+O3IOO+O29OO+O32OO+O3IOO+O32OO+O25OO+042OO+O32OO+O33OO+O2 400N 
+01000+02000+01000+02000+03000+01000+01000+01000-01000-01000-02000N 
+02000=01000+02000-02000+0 3000+04000+02000-02000+06000+01000+02000N 
-02 390-02250+0 32IO-OH2O+OO99O-OO2IO+OI25O+O625O-O235O+O325O-O249OH 
+OO632+OIIH-O2IO3+OI246-OOO2 9+00632+02567-032 i i-8+01112+01113=02468N 

J 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



KSL 5°50 - 26l 
DELETE ROWS AND COLUMNS FROM A MATRIX 

Given a matrix of R (rows), by 6 (columns) size, this 
routine extracts and prints out a submatrix r (rows) 
and c (columns) deleting those rows and/or columns 
designated by the user. The input matrix is a rectangular 
matrix of signed elements as is the output submatrix • 
The size of matrices involved must satisfy the inequality 
(R-M) + (C-c) + C < 822. 

Stops 

1. Master tape KSL 5-50 3^070 

2. Specification tape 2^087 

3. Data tape 2koQ7 
Step 3 may be repeated, if you wish to delete the same 
rows and columns, with a new data tape by raising black 
switch* To read in a new specification tape, insert 
the new tape in the reader and raise the white switch up 
and down. The computer then stops at 2*1-087 for a data tape. 

1. Master - duplicate library copy KSL 5»50 

2. Specification tape - this tape consists of three parts; 

a ) a set of signed integers to indicate the rows R to 
be deleted terminated by an N. (A row is a group of 
C elements terminated by an N). If no rows are to 
be deleted, only an N is punched on tape. 

b) a set of signed integers to indicate the columns C 
to be deleted, terminated by an N. If no rows are 
to be deleted, only an N is punched on tape. 

c) an unsigned number followed by a fifth-hole character 
to indicate the number of decimal places desired in 
the results. 

3. Data tape - the data tape consists of rows (R) of signed 
fractions (c); each row terminated by an N and the 
final N terminated by a J. If an F is used instead 
of an N, the computer stops on a 3^090, allowing another 



tape to be inserted into the reader and read with a 
black switch. This allows the user to break a large 
tape into smaller, easier handled tapes, 
TIME ESTIMATE Time in secondss (.006) RC[d Q + l] + (.05) rc[d Q + 1] 

FF STOPS FF location 0S7 - sum check failure on master readin. 

FF location 092 - number of elements in a row do not 
agree with number in first row. 



- 2 



PURPOSE 



KSL 5 .52 
MATRIX INTERLEAFER (Augment a Matrix) 

This routine takes a matrix A into the machine 3 stores 

it and upon reading in a matrix B punches a new matrix 

C out f Matrix C is composed of row 1 of A + row 1 of 

B terminated by an N 3 row 2 of A + row 2 of B terminated 

"by an ; N| etc.^ until all rows of A + B have been punched 

out. The number of digits in one element of the elements 

of A need not be the same as those in an element of B. 

No terminating symbol need be employed at the end of 

rows A :or B. There must be a terminating symbol at the 

end of each of the two matrices. .The number of columns 

of A or B must not exceed 810. The number of digits in 

one element of A or B must not exceed 12 . Matrix A must 

satisfy this inequality (c + l). r < 102*4-0. To augment 

a a "~* 

a matrix for K-17* K-8., etc., the output will be A -f -B^ 
with an N at end of row. To prepare an augmented matrix 
for M-13, Mr2ky etc., the user can direct that the augmented 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



matrix have an N after row A. and on N after row B^» 



Master KSL 5.52 
Parameter 
Matrix A 



Stops 
2400F 
2te5 
20026 

ko Matrix B 24O0F 

To read in new parameter 3 insert parameter and black 
switch up . Then., read in matrix A and B with black 
switch. At the 2^00F stop a white switch up and down 
will copy the tape from the reader until a fifth-hole 
delay is reached. 

1. Master - duplicate library copy KSL 5«52 

2. Parameter - consists of C space R space d space 

a a a 

c, space d, space x 
where 
G = number of elements (columns) in a row of matrix A «• 



R = number of rows in matrix A (will be same as number 

a x 

of rows in matrix B). 
d = number of digits to be punched in one element of 

matrix A o 
C = number of elements (columns) in a row of matrix B. 
d, = number of digits to be punched in one element of 

matrix Bo 
x = 1 if it is desired to punch an N after each row of 

matrix A and B in augmented form. 
x = if matrix output is desired as an augmented matrix 

with no N at end of rows of matrix A, 
J. Data A - this is a matrix (subject to capacity- 
limitations - see PURPOSE) of signed fractions. Each 
row of the matrix may or may not be terminated by an 
B s J s F or Lo However the last row of the matrix 
must be terminated by an N 5 J, F or L (or + or -). 
4» Data B - prepare same as Data A. 

TIME ESTIMATE ^ k& A R A c + ^(d + l) v B c 

= — ■ — ■ — " ~ — ™ seconc ^ s ( see parameter for 

explanation of 
symbols ) • 

FF STOPS FFOOO - sum check failure on master read in - reread. 

FFOLO - drum transfer failure or input could be incorrectly 

prepared » 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



TIME ESTIMATE 



KSL 5.5^ 
PRODUCE DATA FOR USE IN M-22 

This routine reads in a lower triangular matrix^ signed^ 
scaled Kf 1 with an N terminating symbol at end of matrix . 
It also reads in parameters needed in M-22 and punches 

output to be used directly in M-22 library routine ~ namely 

-2 

on upper triangular signed matrix s scaled 10 and having 

N s s after each row. The matrix is followed by the parameters 
needed for M-22 (supplied by user). See writeup of M-22 
to determine capacity limitations of this program . 

Stops 

Master KSL 5*5^ 3^0NN 

Parameter 24JK 

Data (to end of A^) 240F9 

M-22 parameters 240NN 

Master must be reread to run new problem. 

1. Master - duplicate library copy KSL 5.5^ 

2. Parameter - d space n space 

where d n = number of decimal places in output 

n = size of triangular input matrix (number of rows) 
3» Data - this is a lower triangular correlation matrix t 
signed 5 scaled 10 "" with an N at the end of the matrix. 
Generally this will be K-2^ K-8j or K~9 output^ but it 
may be a matrix punched on tape by the user. 
M-22 parameters - these are the parameters as needed for 
running the output on M-22 (see that writeup for more 
detailed description). It will consist of n numbers of 
0,, 1 9 or 2 plus 2 sexadecimal characters indicating the 
decimal output of eigenvalues and eigenvectors. All of 
this is followed by an L. 
Essentially readin-readout 



Read in time in seconds = 



(n)(n+l) x d 

500 



d„ = number of 
decimal 
places in 
one element 



(n)(n+l) x d Q 

Read out time in seconds = ■ ^r — d = number of 

decimal 
places in 
one element 

FF STOPS FF stops on master read in - sum check failure 

FF001 - No N at end of triangular matrix. 



KSL 5.55 
TRIANGULAR MATRIX TO SQUARE (or Rectangular) 



PURPOSE 



OPERATOR ' S INSTRUCTIONS 



PREPARATION OF TAPES 



This routine takes a matrix in triangular form (K-8_, 
K-17* K-9 output) and punches it out in square form 
with an N terminating each row and a J at the end of the 
matrix. The user may also delete designated rows or 
columns from the triangular matrix before it is punched 
out* The number of decimal places in each element of 
output must be less than or equal to 8, The size of 
the symmetric matrix cannot exceed 111 rows or columns . 

Stops 

1. Master tape KSL 5»55 3^09N 

2. Parameter tape 2^0NF 

3. Data tape 2^09N 

A new parameter may be read in by raising the black 
switch. 

1. Master - duplicate library copy KSL 5»55 

2, Parameter - if a square matrix is desired (no deletions) 
the parameter is d space n space space 

where d * = number of digits in one element (must be < 8) 
n = number of rows or columns in triangular input 
matrix (< 111) 
If a rectangular matrix is desired , (i.e., columns or 
rows are to be deleted) the parameter is d n space n space 
x space s, space s space . . . s space 
where d n = number of digits in one element (< 8) 

n = number of rows or columns in triangular input 

matrix (< 111) 
x = number of rows to be deleted ( < x < 111) 



'!• 



selected row numbers to be deleted 



(NOTE the number of s.'s = x). 

3» Data tape - the data tape consists of a triangular 
matrix in the following order: 

r 



r ll r 21 r 22 r 31 r 32 r 33 



r r 
nl n2 



nn 



Each number is punched as a sign followed by d_ digits. 



TIME ESTIMATE To read matrix and punch out a square matrix 

n 2 (.020d Q + .037) 

To delete rows and punch out a rectangular matrix 

n (.0025d + .0035) + nx(.Ol67d + .0332 ) seconds. 
See parameter for definition of n d n and x) 
FF STOPS FF123 d Q > 8; s ± =0; or s . > n (computer prints out- which 

of three has occurred). To correct read in correct 
parameter tape with white switch up and down. 
FF000 sum check failure on master - reread master. 

* A triangular matrix of unsigned single digits may be 
read in if d~ is set = 0. 



KSL 5*57 
PAGE OUTPUT CORRELATIONS 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



This routine reads in a triangular correlation matrix 
(+1000 ! s in diagonals)., punches these back onto tape, 
unsealed , and in page form, -with the columns and rows 
correctly labeled with column and row numbers running 
from 1 through v where v is the number of variables in 
the correlation. matrix, v must lie between 2 and 1^8. 
The number of digits in one element of the correlation 
matrix read in must lie between 3 and 9* The output 
from this routine cannot be read back into Illiac. 

Stops 

1. Master tape KSL 5-57 2401K 

2. Data tape 2^4-056 

At this point raise the black switch (if complete data 
has been read in). If 2^056 occurs and there is more 
data to be read in,, raise white switch up and down to 
complete reading of data, then raise black switch once 
all data has been read in. After the results have been 
punched , computer will stop on 2^01K. To begin a new 
problem, insert next data tape, and raise black switch,, 
repeating steps described at beginning at this paragraph. 

1. Master - duplicate library copy KSL 5.57 

2, Data tape - the data tape consists of the output from 
routines K-8, K-9<> KSL 2.40 or any other triangular 
matrix of correlations . 

The machine records the number of fifth-hole characters 
after the first diagonal entry (first +1000 above ). 
Whenever the number of fifth-hole characters are following 
diagonal entries (+1000 ! s) exceeds the number after the 
first, the machine will stop on 2^056. (An extra carriage 
return after a +1000 will cause this stop to occur). If 
the complete matrix has not been read in, raising the 
white switch will cause the machine to continue reading 



until the next time an extra fifth-hole character is 
encountered, (in general,, if your data tape is output 
from any of the above routines and is one tape^ and has 
not "been tampered with), the above stop will not occur 
until all the correlations have "been read in. However } 
there will generally be a short part of data left - the 
means and standard deviations. Be sure to note to the 
operator that "all the data will not read in." 

TIME ESTIMATE In seconds .025 n (d.+l) 

FF STOPS FF location 0J7 sum check failure on master readin - try 

again . 
FF location OIF no + for first diagonal +1000 . 
FF location 021 no 1 for first diagonal. 
FF location 0^-9 no + for subsequent diagonal +1000. 
FF location Ok-E no 1 for subsequent diagonal +1000. 
Raising the white switch for any of the last four FF 
stops causes the computer to continue reading;, your output 
will not be affected. 

NOTES In the special case when d. ■= 3 and n > 99 * "the column 

heading will be out of line for column heading > 99 • 
Otherwise j the results will be punched correctly. 



- 2 - 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



TIME ESTIMATE 



FF STOPS 



KSL 5- .60 
MATRIX NORMALIZATION 

This routine will normalize (convert into a unit vector) 
a matrix by rows or by columns. If r is the number of 
rows and c is the number of columns , and you normalize 
by rows , the capacity is 1 < c < 701 and r is not limited. 
When normalizing by columns 1 < c < 351 and r(c+l) < 102^0. 
The elements of the normalized matrix will be printed by 

rows • 

Stops 

1. Master tape KSL ^.60 3^016 

2 . Parameter tape 2001K 

3. Data tape 0F000 

To begin a new problem at OF, move white switch up and 
down. Computer comes to a 2k0l6 stop. Raise black switch 
to read new parameter and new data. 

1. Master - duplicate library copy KSL 5.60 

2. Parameter tape - +d Q N or J. Where +d Q is a signed 
number < 12 indicating the number of decimal places 
desired in one element of the output. If an N is 
used , normalization will be done by rows. If a J 
is used, normalization will be done by columns. 

3. Data tape - the elements of the matrix are punched 
by rows as signed fractions, with an N after each 
row, and a J after the final N. If, instead of an 
N, an F or L is punched at the end of a row, the 
computer will stop. This enables you to read in 
another tape containing the remainder of the matrix 
by raising the black switch. 

Read in of master - 20 seconds 

Read in, calculation and punch out of normalized matrix 

in seconds - .022 rc(d +l) 

FFQ01 wrong number of elements in a column 

FFQ02 wrong number of elements in a row 

FF003 sum check failure on master read in, try again. 



KSL 5.7O 
SUMS OF ROWS AND COLUMNS OF A MATRIX 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



TIME ESTIMATE 



This routine reads in a matrix, punches out the sum of 
rows and columns . Capacity, 2 times number of rows + 
number of columns < 920 * 

Stops 

1. Master tape KSL 5-70 2*K)5J 

2. Parameter tape 2k0kQ 

3. Data tape ^OkQ 

Data tapes of identical size may he read in with black 
switch c To change parameter, white switch parameter in, 
read in data tape with black switch* 

1. Master - duplicate library copy KSL 5. 70 

2. Parameter* - 005+ 

OOcF OOcF c is number of columns 

007+ 

OOrF OOrF r is number of rows 

21+72N 

3. Data - the data tape is punched with N's at end of 

each row. Scale all data at least 10" to keep 

overflow from occurring • No J at end of data. 

r x c x d 
Input time in seconds = 



Output time in seconds = 



150 



rd Q + cd Q 



d. = decimal places 
in input 



50 



d = decimal places 
in output 



FF STOPS 
NOTE* 



None listed. 

The number of decimal places in output will be six (6) 
unless a change is added to beginning of parameter to 
indicate a desire for more than 6 places. This change 
will be in the form 003+ followed by regular 



OOF 00(d Q )F 



parameter o 



KSL 5.71 
SUMS OF SQUARES OF ROWS AM) COLUMNS OF A MATRIX 



PURPOSE 



OPERATOR ' S INSTRUCTIONS 



PREPARATION OF TAPES 



TIME ESTIMATE 



FF STOPS 



This routine forms the sums of squares across rows and 
down columns of a matrix and in addition it calculates 
the total sums of squares (2R+C). Twice the rows plus 
tne number of columns < 920. 

Stops 

1. Master tape KSL 5.71 2^01^ 

2. Parameter tape 2401S 
3 » Data tape 2^01S 

A black switch up will cause new data tapes to be read 
in on same parameter. For a new parameter _, white switch 
up and down .- read in new data with black switch. 

1. Master - duplicate" library copy KSL 5.71 

2. Parameter - d space - where d = number of decimal 
places wished - 6 or above and < 12 is usually wise 
to take care of scaling. 

3. Data - rows of a signed matrix terminated by N's with 

-2 



final N followed by a J. Preferably scaled 10 
1 



In 



most cases 10 
Input time in seconds = 



is O.Ko 

r x c x d. 



Output time in seconds = 



150 



rd Q + cd Q 



50 



FF stops occur for overflow, 

FF stop location 02J - wrong number of elements in a row„ 



K7 - 188 
ITERATIVE ESTIMATION OF COMMUNALITIES : PRINCIPAL AXES METHOD 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



Given a lower triangular correlation matrix^ this routine 
makes a principal axes factoring the matrix,, re-estimates 
the commonalities ^ substitutes the new estimates ft>r the 
old j, and prepares for another cycle » The latent roots and 
new communality estimates are punched after each cycle, 
The size of the correlation matrix is 18 or less« 

Stops 

1. K7 master (part l) 2*K)l4 

2. Data 3^W1 

3« K7 master (part II from where 
removed in l) black 
switch down OF 

If it is desired to interrupt run before OF is reached ^ 

center black switch. When Illiac stops on 2^0F9., white 

switch up and down. The output will be used in preparing 

the data tape for further iterations . 

1. K7 master (duplicate Illiac library copy) 

2. Data tape - this tape consists of the parameter 
followed on the same tape by two copies of the lower 
triangular half of the correlation matrix. The parameter 
is d JvFf L where cL = decimal places desired in output 

(scaled 10 )| v = number of variables or rows in correla- 
tion matrix and f is the number of factors you desire to 
have communalities estimated for. The correlation matrix 
must be scaled so the sum of squares of elements over the 
entire matrix is less than 1/2. The ones in the diagonal 
of the correlation matrix may be replaced by communality 
estimate < one to speed time of calculation. The correla- 
tion matrix is punched by rows (lower left triangle). 
Each element is signed and the last element is followed by 
an No This matrix is duplicated .= that is^ two copies 
are needed on the data %ape. 



TIME ESTIMATE Dependent on number of iterations and size of matrix. 

11 variables - 1 iteration - 1 minute 
18 variables- 1 iteration - 5 minutes 

FF STOPS None listed. 



K8 - 189 

PRODUCT MOMENT CORRELATIONS/ VARIANCE,, 
COVARIANCES,, MEANS AND STANDARD DEVIATIONS 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



Given a raw score matrix of v variables or measurements 
and s samples or observations this routine gives product 
moment correlations and/or variance-covariances and/or 
standard deviations . The maximum number of variables 
allowed is ikk, There is no limit on the number of 
observations. However,, ±f 3 upon consulting the time 
chart enclosed you find your problem will run longer 
than 2 hours it is best to use K-17 with the interrupted 
run option. If the output asked for in this routine 
is correlations in triangular form 5 the output will be 
usable directly in KSL 1.20^ 1,21, 1.^2 as well as 
other routines . 

Stops 
lo K8 master (operator's copy) 3^-OlK 

2 . Parameter 2^0J9 

3. Data 2401K 

A new parameter may be read in with a black switch. This 

is followed by a new data tape. 

lo K8 master (operator has copy - no need to duplicate) 

2, Parameter d Fd Ld JsSx 

c vc n 

where d = decimal places in correlation output (results 



are scaled 10 ), 



d = decimal places in variance-covariances output 
(do not get this if correlations are to be 

read back into Illiac - set d = 0). 

vc ' 

d = decimal places in means - standard deviations 
m 

s = sample size or number of observations « 

x = if triangular output is desired. 

= 1 If square symmetric output is desired. This 
output will be scaled 10 "" and printed by columns , 
each column being terminated by an N. 



TIME ESTIMATE 
FF STOPS 



3. Data - each observation is punched as a signed fraction 
of up to 12 decimal places • An N must "be punched at the 
end of each row. If an F follows a row, the machine will 
stop and another part of the data tape can be placed in 
the reader and the problem continued. No scaling is re- 
quired on the data other than as described below; If your 

data consists of 3 rows of 3 variables +1 +21 »;321 

+16 - 2 -1192 
-111 + 4- 26 

it would be punched as follows: +001 +21 -0321N 

+016 -02 -1992N 
-111 +00 +0026N 

See accompanying chart. 

FF017 - sum check failure on master read in - reread master 

FF015 - number of variables in a row does not agree with 

number in first row - return problem 

FF016 - variance is negative - return problem 

FF023 - arithmetic error _, possible from scaling. 



- 2 - 



K-6HUXBIAM 



The following tisats are estimates only. They are for use as m guide to help prorid* closer estimates in 
the time required to run I£-8 programs through Uliae. In general, estimates b a s e d on theee times should not 
hm more than + 10fj in error. 

The chart belonr gives the time for fixed intervals of variables In eeeonde per sample, I.e., for a 100 
sample sis* problem, multiply the tine in column tvo by 100 to obtain the time In seconds the problem is to 
run before output occurs. The remainder of the chart gives the time, in minutes, that a square or diagonal 



matrix would require to output, for 


convenience sake, these times an 


> given ft 


* 3, k, 


5, and < 


>, decimal 


places of output. 
















RMffi-IH TIME 


PRUT-OUT TIKI IS MJJlfTJS! 
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TRIAH0ULAR MATRIX 




BQiiSK matrix 








3 dec. k dec. 5 dec. 6 dec. 


3 &ec. 


k dec. 


5 dec. 


o dec. 


5 


1 


1 1 1 1 


1 


1 


1 


1 


10 


1 


1111 


1 


1 


1 


1 


15 


2 


1111 


1 


1 


1 


1 


20 


2 


1111 


1 


1 


1 


2 


25 


3 


1 1 1 1 


1 


2 


2 


2 


30 


k 


1 1.5 1.5 2 


2 


3 


3 


k 


35 


5 


1.5 2 2 2.5. 


3 


k 


k 


5 


ko 


6 


2 2.5 2.5 3 


k 


5 


5 


6 


k5 


7 


2.5 3 3 3.5 


5 


6 


6 


7 


50 


7.5 


3 k k k.$ 


6 


8 


8 


9 


55 


9 


* k.i 5 5.5 


■ 

8 


9 


10 


11 


50 


11 


k.5 5.5 6 6.5 


9 


11 


12 


13 


65 


13 


5.5 6.5 7 7.5 


11 


13 


Ik 


15 


70 


1* 


6.5 7.5 8 9 


13 


15 


16 


18 


75 


16 


7.5 8.5 9.5 10.5 


15 


17 


19 


21 


80 


18 ! 


8.5 9.5 10.5 H.5 


17 


19 


21 


23 


85 


20 


9.5 11 H.5 13 


19 


22 


23 


26 


90 


22 


10.5 12 13 1^.5 


21 


24 


25 


29 


95 


2fc.5 


11.5 13.5 15 16 


23 


26 


29 


32 


100 


2T 


13 15 17 18 


25 


29 


33 


35 


105 


31 


Ik 16.5 18 19.5 


27 


32 


35 


38 


110 


3* 


15.5 18 20 21.5 


30 


35 


39 


te 


115 


37 


17 19.5 21.5 23.5 
18.5 21 23.5 25.5 


33 


37 


k2 


k6 


120 


ki 


36 


kl 


k6 


50 


125 


kk 


19.5 23 25.5 27.5 


38 


*5 


50 


5* 


130 


kS 


21 25 27-5 30.5 : 


kl 


*9 


5% 


60 


135 


52 


22.5 26.5 29.5 33 


kk 


52 


58 


65 


1»I0 


56 


25 28.5 31-5 35 


k9 


56 


62 


69 


iR 


a 


27.5 31 3* 38 


5* 


61 


67 


75 (M.P 



KL1 - 192 
QUARTIMAX ORTHOGONAL ROTATION OF FACTORS 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



TIME ESTIMATE 



The routine finds the orthogonal rotation which maximizes 
the variance of the squared factor loadings . This can he 
shown to be equivalent to maximizing the fourth powers of 
the loadings . Two factors are operated one at a time . 
The capacity is F v < 7^5 where f = factors^ v = variables 
with an N at the end of each column. For a more complete 
description of the mathematical method use refer to the 
library routine writeup. 

Stops 

1. KL1 Master 2^3F7 

2 . Parameter 2407N 
3» Data (black switch down) OF 

If continuous , punching has not begun when estimated time 
has elapsed , raise the black switch to OBEY (center). 
When Illiac stops on 2^090 , raise white switch up and 
down. Punching will ensue. 

1. Master - (duplicate Illiac library copy) 

2. Parameter - 006+ 

OOvF OOyF 

0027+ 

OOfF OOfF 

222 78N 
where v is the number of variables and f is the number of 
factors desired. 



3. Data - a matrix tape of f columns of v variables 



Each 



,-1 



Each column of factor 



element is signed and scaled 10 

is terminated by an N. 

»36fv + cf seconds where c and f can be approximately 

determined from the following tables 



No. of factors 


k 


6 


8 


12 


16 


c 


k 


7 


12 


21 


36 


f (in seconds) 


17 


30 


50 


120 


200 



FF STOPS 



(Centroid factors will usually require a larger £ than 
will principal axes factors). 
None listed. 



KLk - 233 
MULTIPLE REGRESSION ANALYSIS WITH TRANSFORMATION 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



Given a score matrix this routine gives first order 
correlations on subsets chosen from the score matrix 
after performing indicated transformations on the scores , 
It also gives a multiple correlation between a dependent 
and independent variables in the subset. Also obtainable 
are covariances, means , standard deviations , standardized 
and unstandardized regression weights and their respective 
standard errors. If it is not desired to receive this 
information^ the proper decimal output parameter may be 
set to 0_, suppressing the output for that particular 
information. The total number of variables on the data 
tape must be < 50* The largest number of variables dealt 
with in any parameter cannot exceed 22. The number of 
observations is unlimited, [if a large number of analyses 
are to be done on one data tape, please use Kl6 - which 
allows the data tape to be stored on the drume KLk re- 
quires the data to be read into the machine after each 
parameter _, KL6 does not] . 



Stops 
2^-012 
200NL 
2^012 



1. Kl6 master (Operator's copy) 

2 . Parameter 

3. Data tape 

A new parameter and data tape may be now read in with a 
black switch for each. 

1. KLk - operator has copy - no need to provide one. 

2. Parameter - consists of four parts. 

sS : = where s is sample size - number of rows in matrix. 

nN = where n is the number of variables in subset. Must 

be 22 or less. 

d F = where d is the number of places to print correlations 

d L = where d is the number of places to print 
cov cov 

covariances 



3--DO+ = where &□„ is the number of places to print 

.Bb -DO 

standardized regression "weights , their standard 

errors and multiple correlations . 

cL J = where d^ is the number of places to print 
Bu Jdu 

unstandardized regression weights , their standard 
errors and standard errors of estimates , 

(Any of the four decimal output numbers may be set to 

zero). 

II Following these six pairs of sexadecimal characters are 

a series of ones and zeros. The ones refer to the position 
of the variables in the subset you wish to hold successively 
independent. For example, in a subset of 20 variables 
from a data tape of 50 variables you wish to hold the 1st, 
7th and 10th variables in the subset (not over the entire 
set of variables) successively indpendent. The series 
of ones and zeros would appear on tape as 1001000001 . No 
delays, spaces, line feeds or other characters with fifth- 
hole punches must appear in this series , The first character 
after the last one in this series must be a delay, space 
or line feed. 

III From this point on until Illiac reads a single fifth-hole 
delay, whatever is on the parameter in the way of identifica- 
tion will be exactly duplicated on the output tape. If 

no identification is desired follow the fifth-hole character 
at the end of part II above with a single-hole delay, 

IV Next follows a sequence of signed integers terminated by 

an N. These integers provide the clue to the transformations, 
if any, wished on the entire group of variables (not the 
subset) in the data tape. Therefore, if you have a data 
tape of 50 variables and are using a subset of 22, there 
will be 50 signed integers ■=- 28 of which will be +0. 
The transformations available are listed below. 
+0 = eliminate variable 

+1 = use variables as Is unchanged 

2 
+2 = variable is squared - x 



3 
+3 = variable is cubed - x 

+k = take square root of variable - x - x > 
+5 = take log variable over 10, :r Q , x > 10 

+6 = take _ arc sin x 

err 
EXAMPLE - You have a 3© variable, 121 sample matrix. 

You wish correlations, multiple correlations and standardized 

regression weights over a subset of 22 variables, holding 

variables (in the subset) 1, 3* and 5 successively 

2 
independent. You wish to transform variable 10 into x 

and 20 into/fx. Your parameter would appear as, 

121S22N5F0L84-0J10101sC0RN 1961° +0+0+l+l+l+0+0+l+l+2-fO4O+l 

P 
a 
c 
e 
+1+1+1+1+1+1h^1^^^4O+1+1^+0+1+1+0-K)+1+140N 

The computer reads in the 36 variables, 121 • row matrix and 
stores it. When it is ready to calculate, it selects 
columns (variables) 3, ' k s §, 8, 9, 10,-13, 1^/15* l6, 17, 
18, 19, 20, 21, 26, 27, 30/31, 3^ 35 to work with, 
squaring variable 10 and taking the square root of variable 
20 before setting up the matrix to work with*. From this 
subset of 22 variables, the first one is held indendent, 
then the third one, then fifth one. 

3» Data tape - the data tape consists of s_ rows of v 
variables. These are signed numbers. Each of the s_ 
rows are terminated by an N. Scaling over variables is 
advisable,, If you have four values for variable one of 
+1 "21+321+6213, column one would appear on the tape as 
+0001 .......... N 

-0021 ......oc. N 

+0321 ...... c .. . N 

+6213 .......... N 

If an F terminates a row instead of an N, the machine 
will stop on a 3012K. This allows another tape containing 
more data to be read in with black switch up. 



TIME ESTIMATE Eead in of master = 35 seconds. , , 

d.sv. v cL_ 

To calculate and printout = [2v - 10 + -Trx — + -=■ — ] seconds 

where v = number of variables in subset, v. = number of 
— — 1 

variables on one row of data tape, d. = number of digits 

in one element on data tape, d^ = decimal places in output 

— 

(one element) and s = number of rows or samples in data 
tape . 
FF STOPS FF011 location 378 - sum check failure on master readin - 

reread master, 
FF012 location 0J4 - row of variables does not agree in 

number with first row - return 

problem . 
FF013 location 06N - drum transfer failure - try again. 
NOTE [If 2 variables are uncorrelated with each other - have 

correlation - routine will hangup on division hangup or 
00000 order] . 



- k - 



K15 - 258 
ANALYSIS OF VARIANCE BY METHOD OF FITTING OF CONSTANTS 

PURPOSE Given a score matrix of v variables and s samples, this 

routine calculates and punches out RAW SUMS OF SQUARES , 
Fitted Constants, the inverse of the cross products 
matrix and "Accounted for Sum of Squares". If it is 
desired., the cross product matrix may be obtained in 
place of the inverse matrix. The total number of variables 
must satisfy the inequalities: 

^~ [3V - K+3] + V < 792 and V < 30 where V = total 
number of variables and K = number of dependent variables . 

There is no limit on sample size unless the accumulated 

39 
sums of squares exceeds 3 • See the Illiac library 

writeup for a more complete description of use of results 

from KL5. 

OPERATOR'S INSTRUCTIONS Stops 

1. KL5 master (Operator's copy) 3^-0^2 
2c Parameter tape 2^084 
3 . Data 240^2 

A new parameter and data tape may now be read in with a 
black switch up for both. 
TAPE PREPARATION 1. KL5 master - operator has copy - no need to provide 

one. 

2 . Parameter - sS vN k + 2 carriage returns identification 
a fifth-hole delay . 

where s = number of rows in data tape - sample size. 

v = number of variables s must be < 30 . 
and k = number of dependent variables . The identifica- 
tion can be any thing for your own convenience 
in identifying the output. If no identification 
is fiesired, punch after the + 2 carriage re- 
turns and a single fifth-hole delay. 
If k is set to v then the inverse matrix will be replaced 
by the matrix of cross-products. However, this cross- 
product matrix will be in reverse order, i.e.., variable 
1 will be last, 2 next to last, etc. 



TIME ESTIMATE 



FF STOPS 



3. Data tape - each observation or sample is represented 
by a row of signed coefficients indicating the relation- 
ship of the observation to the constants to be fitted • 
The dependent variable/variables is the last one to be 
punched in the row. These signed coefficients .describe 
the complete classification. of each observation; use a 
zero if the observation did not occur in a class and a 
1 if it did. Each row is terminated by an N, If weights 
other than 1 are desired for any row of this observation 
matrix, precede the row by the numerical weight j punched 
as +jJ« The machine in effect treats these rows as if 
they occurred j times. If a terminating N is replaced 
by an F, the machine will stop so that more data can be 
read in, 

EXAMPLE - 10 observations k variables 3 subgroups of 
observations = 3^ 3> and k» 

Parameter tape - 10S7N1 + (2 carriage returns) IIA* 
Data tape - +1+0+0+0392 -0^26+1123-2 3N 

+^J+l+0+0-0263+1023-002^+19N [Treat as though four observa- 
+1+0+0+12 35-0032+1926+03N 
+0+1+0+0032+0092 -0326+09N 
+0+1+0-0069-1023+0.^62 -ION 
+0+1+0+0132+62^5+1111-11^ 

+2J+O+I+O-O294-O236+O232-62N [Treat as though two observa- 

tions 
+0+0+1-2 369+0326+0436-03N 

+0+0+1+2469-0012+0626+05N 

+0+0+1-0235+002*H-0727-03N • 

Read in 10s (v + v) 10" 5 

Compute 2[v-kr 

Output 136 [v 2 + 3v + k - k] 

FF010 location 03K = drum transfer failure try again 

FF201 location 099 = accumulated number exceeds 39 - 

return problem 



2 - 



FE200 location 08F - number of variables in a row does 

not agree with number in first rovr« 
Return problem <> 
NOTE If k > I* then the covariance model is in effect and the 

constants printed out for the k variables are in reverse 

order! i«e«j the last variable first, next to the last 

2nd 5 etco 



- 3 - 



kl6 -.263 
modified multiple regression analysis 



PURPOSE 



OPERATOR ' S INSTRUCTIONS 



TAPE PREPARATION 



TIME ESTIMATE 



FF STOPS 



See "writeup of KL4 for complete instructions on running 
this problem? The only difference is in (l) the data tape 
need only be read in after parameter 1 and (2) the parameters 
have one additional terminating symbol (see below) » In 
addition the sample size times [the total number of variables 
on tape +1] must be less than 8137° 

Stops 
lo Kl6 master (operator's copy) 24-012 
2o Parameter (type N) 24113 

3. Data 24012 

4. Parameter (Type d) 24012 

Any number of parameter (Type' J ) may be read in with a 
black switcho These all will -work on the data read in in 
step 3° To read in a new data tape, black switch in a 
parameter (Type N) and then read in the new data tape., 
This may then be followed by parameters (Type j). 
lo KL6 master (Operator has copy - no need to provide 

one)« 
2. Parameter (Type N) - identical to one described in 

KL4 except an additional N is added after the final N- 
3° Data - as described in KL4„ 
4, Parameter (Type J) - identical to one described in Kl4 

except one additional J is added after the final N. 
[NOTEs transformations for KI6 are only +0, +1, +2, +3, 
+4, and +5« No +6] o 
Read in of master = 35 seconds 

To calculate and print out = [6v + v d ] seconds where 
v = number of variables in subset and d n = number of 
places in output e 
FF011 location 378 - sum check failure on master readin - 

reread master o 
FF012 location 12J - number of variables does not agree 

with number in first row - return 
problem . 



FF013 location 06N - drum transfer failure - try again. 
NOTE If two variables are uncorrelated - have zero correlation 

with each other - machine will hangup on a division 
hangup or a 00000 order. 

See KL4 EXAMPLE. The only change would he an additional 
N at end of Parameter A - an additional J at end of 
Parameter Bo 



2 - 



KL7 - 509 

PRODUCT MOMENT CORRELATIONS, VARIANCES GOVARIANCES, MEANS 
AND STANDARD DEVIATIONS 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



This routine takes a raw score matrix, row by row 
providing center correlations between the columns in 
the matrix. It differs from K-8 in that a problem can be 
completed in two or more distinct machine runs, (See 
"INTERRUPTION OF PROBLEM"),, The number of variables 
must be < iMw The number of rows or sample size is 
unlimited a The output, if obtained in triangular form, 
is ready to be used directly as input f or KSL 1,20, 
.1.21, .1,52. 

Stops 

1. Master KL7 2^01K 

2. Parameter 240J9 
3o Data ZkOUH 

A new parameter may be read in with a black switch. A new 
data is then read in with a black switch. If it is desired 
to run the problem with more than one machine run, the 
following is used as OPERATOR ' S INSTRUCTIONS . 



FIRST RUN 


Stops 
2401K 




1. Master KL7 




2 . Parameter 


240J9 




3« Data tape, part 1 (2N's at end 

of tape) 


FF015 




h. Move white switch up and down 


OF 




(Do not print output) 






SUBSEQUENT RIMS 


Stops 




1. Master KL7 


2^01K 




2. Tape (Output from 1st run) 


240*1-7 




3. Data tape, next part 


2401K 


(If total 



problem is finished. If there is still a third data 
tape, instructions for steps 3 and k are as in FIRST 
RUN instructions). 



TAPE PREPARATION 



TIME ESTIMATE 
FF STOPS 



1, K17 master - (duplicate library copy) 

2. Parameter - dJF d _ L d JsSx 

cO m 

where d~ - = number of places in correlations « 


d n = number of places in covariances. 



m 



number of places in means and standard 
deviations (the above results are scaled 



10 



-1 



if 3 significant places desired set 



d = k). 
c 

s = sample size or number of rows, 
x = if symmetric matrix of correlations and/ or 
covariances is to be punched in triangular 
form with an N at end of the matrix. Use 
if you desire to use output in KSL 1.20, 
1.21, 1.52- 
= 1 if a square matrix is desired. If 1 is used 
correlations and/or covariances are punched 
out by columns with an N at the end of each 
column. 
3. Data tape - the data tape consists of s rows of n 
elements each. The n variables are punched as signed 
fractions . Each row is terminated by an N . If an F 
terminates a row instead of an N, the computer will stop 
on 300^9. To continue reading the remainder of the data 
tape, raise the black switch, being sure the new tape 
is in the reader before doing so. 

If you wish to run the problem in interrupted machine 
runs (see OPERATOR'S INSTRUCTIONS), at the end of the 
first part, punch an extra N terminating symbol. 
See accompanying chart. 
FF001 location 152 - sum check on reading sexadecimal 

tape has failed. Reread master, try 
sexadecimal tape again. 
FF015 location O^N - an extra N has been punched following 

a normal N. White switch up and down 
to get sexadecimal tape for next 
run .- or the number of variables in 



_ 2 - 



one row does not agree with the 
number in the first row. Raise 
the white switch up and down to 
get partial correlations to this 
point - correct row in error - run 
rest of data as instructed for an 
interrupted run* 
FF016 location 088 - the variance is negative and outside 

the tolerance limits. Remove problem « 
FF017 location ^¥h - sum check on reading master tape 

has failed . Reread master, 
FF025 location 08l - arithmetic error possibly due to 

scaling on data tape. 
DATA OUTPUT If both correlations and covariances are specified by the 

parameter tape., the correlations appear in the first column 
and the covariances in the second column. The correlations 
are scaled 10° • The scaling of the covariances depends 
upon the scaling of the variables on the data tape. If x. 

is scaled by 10" and x. is scaled by 10 , then the 

J - (i+i ) 
covariance i j is scaled by 10 . 



time mnwBp% jor k-6 program 



The fallowing times are estimates only. They are for us® as a guide to help prowl** closer estimates in 

tiro required to run K-6 progress* through Illiac. In general, estimates based on these times should not 

be mare than ♦ 10)1 in error. 

The dart below gives the tims for fim»d intervals of variables In seconds per sample, i.e., for a 100 
MBftjle sice problem, multiply the tims in column tuo by 100 to obtain the time in seconds the problem Is to 
nmbefore output occurs. The remainder of the chart gives the time, in minutes, that a square or diagonal 



itrlx liould require to output, 
places of output. 

BEAD-IS TINE 
10. OF VARIABLES fflBD/SAMK£ 



Tor convenience sake 



these times are given for 3, K 5, end 6, decimal 



5 

10 

15 
20 

25 



30 

35 
40 

45 
50 



55 
50 
65 
70 
75 



80 
85 
90 
95 

100 



105 
110 

115 
120 

125 



130 
135 

l40 

1jS4 



l 
l 

2 
2 
3 



4 

5 

6 

7 
7.5 



9 

11 

13 
14 
16 



18 

20 

22 

24.5 

27 



31 
3* 

37 

4l 

44 



48 
52 
56 

a 



PKEHT-00T TIME IH MDTOTHB 
TRIAHOULAR MATRIX 
3 dec. 4 dec. 5 dee. 6 dec. 



1 

1 
1 
1 

1 



1 

1.5 
2 

2.5 
3 



1.5 
2 

2.5 

3 
4 



4 

*.5 

5.5 

6.5 

7.5 



k.5 
5.5 

6.5 
7.5 
8.5 



8.5 

9.5 

10.5 

11.5 

13 



9.5 

11 
12 
13.5 
15 



14 
15.5 

18.5 
19.5 



16.5 
18 

19.5 

21 

23 



21 

22.5 
25 
27.5 



25 
26.5 
28*5 
31 



1 
1 
1 
1 
1 



1.5 
2 

2.5 
3 

4 



5 
6 

7 
8 

9.5 



10.5 

11.5 

13 

15 

17 



18 

20 

21.5 

23.5 

25.5 



27.5 
29.5 
31*5 
3» 



1 
1 
1 

l 
1 



2 

2.5 
3 
3.5 

*.5 



5.5 
6.5 
7.5 
9 
10.5 



n.5 

13 

14.5 

16 

18 



19.5 
21.5 
23.5 
25.5 
27-5 



3 dec. 



2 

3 

4 

5 
6 



8 

9 

11 

13 
15 



30.5 
» 
35 
38 



17 
19 
21 

23 
25 



27 
30 

33 
? 

38 



41 
44 
*9 
5* 



8<JJARX MATRIX 

4 dee. 5 dec. 6 dec. 



1 
1 
1 
1 
2 



3 

4 

5 
6 
8 



9 

11 

13 
15 
17 



19 
22 
24 
26 

29 



32 

35 
37 
41 

%5 



52 
56 
61 



1 
1 
1 
1 
2 



3 

4 

5 
6 
8 



10 
12 
14 
16 
19 



21 
23 
25 
29 
33 



35 
39 
42 
46 
50 



5* 
58 
62 



1 
1 
1 
2 
2 



4 

5 
6 

7 
9 



11 
13 
15 
18 
21 



23 

26 

29 
32 

35 



38 
42 
46 
50 
5* 



60 
65 
€9 
75 



C*. 



L7 - 230 
AUTOMATIC LINEAR EQUATION SOLVER 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



The routine takes a set of simultaneous linear equation 

coefficients set up in matrix form (A and solves the 

* nn 



unknowns x n , x» , 



X- 



x 



If desired it will also 



3 n-1 
calculate and provide the residuals » The values for the 

unknowns are punched out in a single column with a decimal 

point properly placed in each result,, The capacity is 

1^3 equations or less. 

Stops 

1. L7 master (operator's copy) 3^-OlN 

2 . Data 2^009 

If residuals are wished at this point instruct the 
operator to re- insert data into reader and white switch 
up and down. If no residuals are wished , read a new 
data tape in with a "black switch up, 
lo L7 master - operator has copy - he will supply it. 

If request is made on OPERATOR'S INSTRUCTION sheet. 
2. Data tape - the elements in the data tape are signed 
fractions 5 scaled so that each coefficient or element 
is < l/2. Each element must have 12 or less digits « 
■ Each row of coefficients must be terminated by an N. 
The number of rows in the data matrix will be one less 
than the number of coefficients. At the end of the data 
matrix,, after the final N 5 place a number d from 2 thru 
N to indicate the number of digits in the elements to be 
punched out. K = 10 places^ S = 11 ^ N = 12. 
An example follows 2 
56 x, + 23 x n + 9 x = 10 

L d j 



■2 x ± + 3 



2 x 



3 



+31 x n 



32 x 2 + 8 x^ 



-5 
, x ^ ^ . „ ^ - 999 

The data tape asking that solution contains 10 decimal 

places in each result would be as follows s 

+O56 +023 4009 -010N 

=2 +3 =2 +5N 

+0031 -0032 +0008 -0999NK 



TIME ESTIMATE a) 15 seconds to input program 

b) .OO^K seconds to input coefficient where IC = (n)(n-l) 

c) .OOlfn seconds to solve equations 

d) .017(d +3)n seconds to punch results 
where n - number of elements in one row, 

FF STOPS FF02K - master sum check failure (read master again) 

FF026 - drum failure during triangularization (try again) 
FF027 - drum failure during black substitution (try again) 
FF028 - wrong number of coefficients (compared to first row) 
FF029 - trying to solve for more than 1^3 equations. 



_ 2 - 



PURPOSE 



OPERATOR ' S INSTRUCTIONS 



TAPE PREPARATION 



AUXILIARY M7 - 150 
PRINCIPAL AXES FACTOR ANALYSIS 

This routine takes a lower triangular matrix of correla- 
tions,, scaled 10*" 2 and prints out up to 23 factors obtained 
by scaling the eigenvectors by the square root of the 
corresponding eigenvalue (latent root). Each factor is 
printed as a column terminated by an N. The eigenvalues 
are punched as a separate column after the last factor, 

Stops 
lo M7 master (part i) 2^00K 

2. Data 2^7 

3. M7 master (Part II - begin where 

removed in l) OF 

The master must be reread for a new problem * 

1. M7 master - (duplicate Illiac Auxiliary library copy) 

2, Data - the first part of the tape consists of the 

-2 
correlation matrix to be factored. It must be scaled 10 

signed and the entire matrix terminated by an N. The 
correlation matrix will be the lower left triangular half 
punched by rows. The second part of the tape consists 
of two numbers each followed by a sexadecimal character - 
they are d JvF where d Q is the number of decimal places 
desired in the output and v is the order or number of 
rows in the correlation matrix. 
An example of data input follows z 
+1000 ) 

+0^38 +1000 ) This could be output from K2^ kS> -K9* 

+0072 =0309 +1000 ) ^ 

The data tape would be 

+01000 

+00^38 

+01000 

+00072 

-,00309 

+01000N 



5J results will have a sign,, zero and four significant 
figures (i.e., scaled 10 ~ ) 

3F order of matrix (number of rows) is three. 

THIS MUST BE ON ONE TAPE 
TIME ESTIMATE (Same as for M3-117) 

FF STOPS None listed. 
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ML2 - 173 
MATRIX TRIANGULARIZATION AND COMPUTATION OF ITS DETERMINANT 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



PREPARATION OF TAPES 



DATA OUTPUT 



SCALING 



Given any number of square matrices up to kO x k-0 elements 
in size, this routine will upper triangularize and compute 
the determinant value - punching either or both outo 

Stops 

1. M12 master 20095 

2 . Data 20095 
Additional data tapes may be read with a black switch, 

1. M12 master - duplicate library copy 

2. Data tape - a square matrix, punched row by row, each 
element signed and scaled so its value is less than 1/2 
(see Scaling )« Each row of elements must be terminated 

by an No If it is desired to print out triangular matrix 

end last row of matrix with a J instead of an ,N# 

To interpret from results the value of the determinant 

use the following: 

(-1)% x 10 2 , ; -1 < b < 1 

Example - 013 -l6l2 3^56789 represents 

if n even -,16123456789 x 10" 15 

n odd .16123^56789 x 10~ 15 

If an upper triangular matrix is to be punched out, it 

will be punched out by columns with eleven digits , 

All elements should be less than l/2. In order to determine 

value multiplied by determinant result to get correct 

results follow example below; 

if original matrix is 1 2 3 

5^1 data tape would 

. .01 o004 .0091 

be +.1 +.2 +.3N x 10 multiply value of 



+.05 +.04 +0OIN x 10' 
+.1 + .04 +.091N x 10 



-1 



not 

on 

tape 



determinant obtained by 

2 

10 to get correct result 



10 



2 -5 

TIME ESTIMA.TE a) Read in matrix [ (d +l)n + n] 10 seconds 

"b) Computing time and punching determinant 
= r n 10 seconds 

c) Punch time of triangularized matrix 

2 -^ 
336 (n + n) 10 seconds 

FF STOPS Only in library subroutines Q 
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ML3 - 179 
COMPLETE LINEAR MATRIX EQUATION SOLVER AND GENERAL MATRIX INVERSION (N < 22) 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



This routine solves the linear matrix equation Ax = B 
for matrices A and B -when A is non-singular and a square 
matrix of order n x n and B is of size n x m. The size 
of n and m must be such that the inequality 
nm + g (n + l) + (n + m) < 842 is satisfied. 
The solution x of size n x m is punched out by columns 
(an N at the end of each column) The last element in 
each column is a scaling factor (see example). If a 
matrix A is to be inverted , the capacity is a matrix 
A of size 22 x 22. 

Stops 

1. Master ML 3 2kO~5F 

2. Data 2403F 

A new data tape may now be read in with a black switch. 

1. Master Ml 3 - duplicate library copy 

2. Data - the data tape for solving Ax = B is prepared 
as follows; Punch row 1 (n) of matrix A followed by row 
1 (N) of matrix Bo Then row 2 (n) of matrix A followed • 
by row 2 (N) of matrix B. Continue until row n (n) of 
matrix A and row n (n) of matrix B are on tape. Each 

row of A and each row of B is followed by an N. The final 
N on the tape is followed by a sexadecimal character 2 
thru S indicating the number of digits desired in the 
output. All rows must be scaled so that no element is 
> 1/2. Each element must be preceded by a sign and 
must have digits < 12 in number. 

The data tape for inversion of matrix A consists of n 
elements and n rows^ each element signed and < 12 in 
number of digits. The elements in the matrix must be 
scaled so no element is > l/2. Each row must be 
terminated by a J. Following the nth J will be a sexa- 
decimal character 2 thru S indicating the number of digits 
desired in the output (inverted matrix). 



TIME ESTIMATE 



FF STOPS 



EXAMPLE (input) 



EXAMPLE (Output) 



Read in a matrix A and B = n ., ,-~ d. seconds 

150 1 

Computing time for solving Ax = B when B is in columns 



[n 3 ] 
500 



m seconds 



2 5 

Computing time: fir inversion of A = kn - «5n + .Olir seconds 



n d 



Punching of output = 







50 



seconds 



No FF stops. However, if an F is punched on the output 
tape with no other output, it indicates either overflow 
of elements or a singular matrix * 
An example would he 



A = 1 


2 





3 


-4 


1 


=2 


5 


6 


B = 


1 


2 







-1 




-2 






The data tape for solving Ax = B would be 

+1 +2 +0N +i:+2N 
+3 -k +1N +0 -IN 

-02-+05 +06n-02 +00NK (last row scaled so that elements 

= 1/2) 

Output will he to 10 places rounded off. 

The data tape for inversion of A would he 

+01 +02 +00J 
+03 -Ok +01J 
-02 +05 +06 JK (all elements scaled alike so no element 

> 1/2) 

A tripical output for a 3 x 3 matrix may look like the 

following; 

+1*000 
-0.500 
+2.500 
+1.000N -» scaling factor - decimal point after 1 

-20.00 
+00.30 
-25 .00 
+01.00N 

+030.0 
+010.0 
+005.0 
+001. ON 



M 21 - 240 
MATRIX MULTIPLICATION 



PURPOSE 



OPERATOR'S INSTRUCTIONS 



TAPE PREPARATION 



The routine allows two matrices A. . and B., to be 

ij ik 

multiplied . j < k-12 and jk < 102^0 . Each element of 
the matrix must "be less than one. If overflow occurs, 
the program records it and punches out the proper 
element of the product matrix C accordingly . This 

1.K 

would be punched as signed integer,, a decimal point 

and an n place fraction otherwise (i.e.., elements < l) 

the elements of the product matrix C, is punched as an 

n place fraction. The elements of C are preceded 

lis: 

by three constants is the number of rows of C., ; k, 

the number of columns , and n the number of decimal 

places printed. Each row is terminated by an N. The 

matrix C can be read back into the computer if no 

elements are greater than one. If you wish to take C , 

ik 

output and place it directly back into the computer in 
the same production problem, it is wise to be sure that 
all elements of A. . and/or B are originally scaled 

n IJ Jk 

10 ' , that is, each element has a leading zero. 



Master M21 



Stops 
30026 



B matrix (parameter on beginning 

of tape) 30030 



3° A matrix (parameter on beginning) 30019 
A new A matrix may be read in with a black switch. If 
you wish to read in a new B matrix and A matrix, white 
switch then black switch in the new B matrix and A matrix, 

1. Master M21 - duplicate library copy 

2 . B matrix - this matrix is preceded by a parameter 
on the same tape. This parameter is j space k space d n 
space, where j is the number of rows of matrix B| k is ■ 
the number of columns and n the number of decimal digits 



to be punched in each element of C 



ik 



The elements of 



B are punched by rows, the elements being signed fractions 



Each row is terminated by an N (replacing an N by a J 
causes the computer to stop., allowing a new tape with 
the remainder of the matrix to be inserted in the reader). 
3° A matrix - punched as B matrix, except the parameters 
are i space j space d j i is the number of rows in A, j 
is the number of columns and d is the same as in the B 
matrix parameter * 

TIME ESTIMATE Input of B matrix if d. is the number of digits in one 

element in seconds (kjk d. + l8jk + 80k + 500 ) 10"^ 
Input of A matrix,, calculation and output in seconds 
(30ij d. + ki$ &. + 6o6ik d Q + 200ik) 10~ 5 

FF STOPS FF03N - sum check failure on master readino 

FF03J - drum transfer failure . (This will occur if the 
j (number of columns ) in matrix A does not agree 
with the j (number of rows) in matrix B) . 

NOTE For instructions on performing a series of A • B multipli- 

cations with just new A matrices or with both new A ma trice 
and new B matrices without stopping see the the more 
complete library writeup. 
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M2^ - 260 

COMPLETE LINEAR MATRIX EQUATION SOLVER AND 
GENERAL MATRIX INVERSION ROUTINE 



PURPOSE 



OPERATOR' S INSTRUCTIONS 



TAPE PREPARATION 



Given two matrices A. . and B. . this routine will solve 

11 ij 

the matrix equation A. . x = B. . for x where A is non- 
11 ij 

singular (its determinant value f 0). It also will take 

a matrix A., and give its inverse A.." • The limits of 

11 P 11 

B. . are set by i + j < l6k and 1 + 2ij + 3i < 20, teO. 

-*- j 
To get the inverse i < 82 

Stops 

1. M24 master 2ko6K 

2 c Data tape 2^06K 

A new data tape may be read in with a black switch up* 

1. M24 master - duplicate library copy 

2o Data tape - when matrix A and/or B are to be read 

from tape 3 they appear as a matrix in this form 



± A 11 


±*12 


-f A 

-13 


000 


+ A n . C 

— li 




± B 11 


± B 12 


000 


±B U C 


±*21 


- A 22 


-•- A 

- 23 


000 


± A 2i C 




± B 21 


± B 22 


000 


± B 2i C 


c « 


o 


» • • 








+ A.. 
— il 


±A 12 


± A i3 


OOP 


■+ A.. C 

— 11 




+ B. -, 
— il 


± B i2 


000 


+ B. . C P 



where J - A . and + B. . are the elements of the matrices o 

- 11 — ij 
C is a control character for calculations. Pis a control 

character for output • 

If the problem is to solve a linear matrix equation,, then 

the C after each row is an No If the problem is the 

inversion of A... C is a J and the elements B. . are not 
11 iJ 

punched . P determines the form of output . If 1 < P < S 
output will be on tape by columns where p is the number 
of digits to be punched for each element . < Each column 
will be followed by a scaling factor or an N (be sure to 
remove the scaling factor before reading the output from 
M2^ back into the computer)* 



If P = L the results of the computation will be placed 
consecutively on the drum "by row, starting at location 
2560. If the user wishes to utilize the drum for series 
of computations, he is referred to the more complete 
writeup of this in the regular 111 lac library routines 
book 
TIME ESTIMATE Inversion: i = 20 = 90 seconds 

± = kO = 600 seconds 
Linear equation: i = 10, j = 1 ~ 8 seconds 

i = 10, j = k ~ 15 seconds 
For higher orders, the time increases approximately as n . 
FF STOPS FF010 - drum transfer failure on y - 1 

FF011 - sum check error in reading program master 
Read again. 
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M28 - 270 
MATRIX MULTIPLICATION WITH OR WITHOUT RESCALING 



PURPOSE 



The routine will form the product of two matrices s 
(A. . B. .). The product ma 
"tape or stored on the drum 



(A. . B. .). The product matrix (C. .) can be punched on 

Matrix B. . can be rescaled 
by any power of ten prior to multiplication A set 



of products of various sizes can be produced without 
rereading parameters s and these products can be stored 
on the drum for use in subsequent operations. Multiple 
products such as (A) p (a)^~ o . . AB = C can be formed 
without punching intermediate results on tape. Capacity 
limitations are found under the heading "STORAGE OF B 
AND C MATRICES ." 
OPERATOR'S INSTRUCTIONS k cases are cited. The first is for a simple (A. . B. . = 



C.) 



multiplication . 




CASE I 


Stops 


1. M28 (Operator's copy) 


3^02L 


2. Parameter tape, type T 


2k060 


3» Data B 


240J7 


ko Data A 


2402L 



TAPE PREPARATION 



To begin a new problem at 2402L read in a new parameter 

with a black switch . 

CASE I 

1. M28 master - operator has a copy at Illiac reader 

console,, Request as above o No master needed if 

requested as above o 
2c Parameter tape, type T - this tape consists of three 

elements and is the same for CASES I. II and III. 



They are d n space + s and x where 
d = decimal output in matrix C . . when d n < 13 



If d, 



is > 17), d n will be interpreted as a drum address 
at which to store matrix C. . 



^ ■" "0 
a drum address 

See section called 



STORAGE OF MATRICES B AND C. 
+s = scaling exponent. Omit this if no rescaling is 
desired on B. .. If +ls is given each element B. . 
will be multiplied by 10 before the matrix is stored 
on the drum. 



If --Is is given, each element of B. . will be multiplied 
" 1 j 

by ,1. If your B. . matrix has a leading zero before 

each element, (i.e., +092 -032). a + Is will eliminate 

the zero. However, if there is a chance that the 

product matrix C. . will have elements over +1.00, 

-*- J 
rescaling matrix B. . would be reverse. A scaler 

of -12 would cause the elements in matrix B to 

all become zero - a thing to avoid. 

x = a directive, either N, J, F or L. When x = J or L 

and d is < 13, the rows of C will be punched in 
u 1J 

columnar form with a carriage return after each 

element. An N will be punched after each row and a 

J at the end of the matrix. When x = N or F and 

d is less than 13, the rows of the matrix C.. will 
iJ 

be printed across the teletype page. The rows will 

be evenly punched so that a neat printout will result. 

An N will follow each row of C. ., and a J will be 

—J 

punched at the end of the matrix. 

Whenever d n is 13 or greater, matrix C. . will be stored 

on the drum, and no distinction is made between the 

directives. For example, a parameter of 10 space n 

means print C. . with 10 decimal places by rows across 

the teletype page<= 

DATA TAPES A . and B. . The elements of both tapes (all tapes in later cases) are 
ij iJ 

punched by rows as signed fractions with an N at the 

end of each row. If any row is terminated by an F 

instead of an N, the computer will stop. By raising 

the black switch, a second part of the matrix will be 

read. A terminating symbol is required at the end of 

each matrix. Matrix B. . can have a J or L. If matrix 

A. . has a J terminating symbol - the next tape read 

iJ 
in (if any) must be a parameter tape. If an L terminates 

an A. . tape, the previous parameters is retained and a 

new B. . may be read in on this parameter. 
U 



OPERATOR'S INSTRUCTIONS (CASE II ) (A p ) (B P ) = C P 

1. M28 master (operator's copy) 
2 



3. 

k. 

5- 
6o 



Parameter tape, type T 

Data tape B^ 1 ' 

Data tape A^ ' (with L terminat- 
ing symbol) 

Data tape B 



Data tape A 



(2) 
(2) 



Stops 
3^02L 

2ko6o 

2taT7 

2ko6o 

240J7 



(with L terminat- 
ing symbol) 



2ko6o 



TAPE PREPARATION 



etc. 
(CASE II) 

1. Master (see CASE i) 

2. Parameter tape., type T (see CASE i) 

3. Data tapes (see CASE i) 



That only difference 



^s that data tapes A^ p ' must 



contain an L as a terminating symbol. 
OPERATOR'S INSTRUCTIONS (CASE III) A^ A^" 1 ' . . . A^ A^"' B = C 

1. M28 master (Operator's copy) 

2. Parameter tape, type T 



3« Data tape B 

k. Data tape A^ ' (j terminating symbol) 

5" Parameter tape, type D 

(2) 

6a. Data tape A^ ' (L terminating symbol) 

6b. Data tape A^' (L terminating symbol) 



TAPE PREPARATION 



Stops 
3^02L 

2ko6o 

2^fOJ7 
2402L 
21+OSJ 
240SJ 
240SJ 
Continue reading Data tapes A^ ' ... A^ p ° I. A^ p ' with 
a black switch. At the end of A^' a J terminating 
symbol will cause the product matrix C to be printed. 
(CASE III) 

1. Master (see CASE i) 

2. Parameter tape - type T, since the first product 
A* ' B must be stored on the drum at a location, 
the d directive will be greater than 13. It can 
be determined by referring to the section below 
titled STORAGE OF MATRICES B AND C. An N or J 



must terminate this parameter. For example, a 



parameter of "^558 space +2N" would cause A 



a) 



- 3 



B to be stored by rows in successive locations at 

drum location k^Q after having multiplied each element 

(1) 

of matrix B by 100 before forming the A K ' B product. 

The type D parameter tape must be terminated by either an 
F or an L. (When this is read in and the black switch 
lifted the product matrix A^ ' B is read from the drum to 

form A^ ' A^ ' B is A^ 'is read into the computer. This 

(1) 

product matrix A x B will be read in from drum location 

D, rescaled if so indicated^ and restored as a new B 
matrix at location 25S0 on the drum). For instance a 
type D parameter might be "6 space -IF". This means read 
rows of previously formed product matrix from the drum_, 

multiply each element by l/lO and store as a B matrix: 

(2) , 
form the new product matrix of A v (this matrix is read 

(l) 

in at this time) A^ ; B and print it out to 6 places as 

rows across the page. 

DATA TAPES (See CASE i) Data tape B and a! 1 ' will have a J terminating 

symbol at end of matrix • All other data tapes will have 

an L terminating symbol until the final product is wished 

printed o The last data tape read in will then have a J 

terminating symbol. 

OPERATOR'S INSTRUCTIONS (CASE IV ) A. M = C. where M is at drum location D 

' lr re ic re 

and Williams Memory has been cleared. 

Stops 
1. M28 master (Operator's copy) 3^02L 

2a. Parameter tape 5 type D 240SJ 

2b. Parameter B - (white switch up and down) 2^+OSF 
3« Data tape A, (with J terminating symbol) 2^02L 
This case would be used in the instance where a product 

matrix M had been formed and stored on the drum - and 

re 

was still present on the drum (i.e., the drum had not 
been clear ) e It would be possible to read in a Data 
tape A. and print out the new product C. . 
TAPS PREPARATION 1. Master (see CASE l) 

2a. Parameter type d - see CASE III 



- k 



2b, Parameter B-- Consists of D space r space c space 
where D is the drum location of M (see STOEAGE OF 



re 



B AND C MATRICES) r is the number of rows in the 

stored matrix M and c is the number of columns 
re 

in M . 
re 

STORAGE OF B AM) C MATRICES 

The first element in each row of B is stored at 2560+ 
(j-l)l where j is the number of rows in matrix B and I 
is determined by the number of columns in B. The capacity 
for B is shown in the chart below: 



Columns in B 


I 


max. rows 


< k < 65 


65 


157 


65 < k < 129 


129 


79 


129 < k < 193 


193 


53 


193 < k < 257 


257 


39 


257 < k < 311 


321 


31 



If matrix C is not stored on the drum^ the number of 
rows of matrix A is not limited. The maximum number of 
columns in matrix B is 311 « 

The first available location for the storage of matrix 
C on the drum will be after the final element of matrix 
B. This can be determined by the following formulas 
D = 2560 + (j-l) I+k (see above paragraph for explanation 
of jy I and k). The number of locations required for 
storage of matrix C will be i (k+l) and D+i(k+l) < 12 .,800. 
For CASE III^ when multiple products are formed D must 
be chosen large enough to accomodate the largest B 
matrix o -As an example: you wish to form a product matrix 



of the following size matrices 



A<3> = 



12 rows 30 columns 1 



A^ ' =30 rows by ^0 columns | A^ ' = kO rows by 15 columns 

and B = 15 rows by 90 columns. This will give a product 

matrix C of 12 rows by 120 columns. The sizes of the 

matrices stored ares 

A^ 'B = J+0 rows x 90 columns = 36OO 

A (2)|- A (l; B j = 5Q rowg x £ columns = 2900 



5 - 



A (3]j- A ( ^Cl; B j _ 22. rows x 90 columns print out* 

A^ B provides largest matrix to "be stored. Even though 

B read in must "be stored, its size (15 x 90 ) is smaller 

(l) 

than A^ 'B to he later stored from drum location 2560 on« 

So matrix A^ ^B will he termining matrix as to size of 
drum area to he used. In this product matrix j = 90 and 
k = kO. 

From the chart if k = k-0, I = 65 and maximum j allowed 
is 157 • Therefore the drum location for storing C matrix 
will he 2560 + (90-1) 65 or 8^5 = D- 

As A^ 'B will also he the largest G matrix stored on the 
drum it will require i (k+l) locations for storage. Since 
i = ^0 and k = 90 this becomes ^>6kQ locations . The 
number of locations for storage of B and C must fulfill the 
inequality D + i(k+l) < 12800. 

83^5 + 36kO = II985 < 12800. Therefore in the parameter 
needing a D drum location D = 83^5 and enough space is 
left for storage of matrices C as they are formed. 
TIME ESTIMATE To read master and parameters - 30 seconds 

To read matrix B from tape and store on drum ° 

jk(o0Q5d. + .013) seconds 
To read matrix C from drum and restore as B matrix - 

i(o005k + .018) seconds 
To read matrix A and calculate - 

ij(.005d. + o010) + ,ik(.00i|-j + *009) seconds 
To punch matrix C - 

ik(.0l67d + .03^) + «067i seconds 
To store matrix C on drum = 

i(*002k + .008) seconds 
FF STOPS FF001 location 06K J or L at end of 1st row of B; white 

switch to continue. 
FF002 location 08^ Number of elements of B are not equal. 
FF003 location 0F8 Row vector (size does not conform to 

column vector of B. 



FFOO^J- location 0S2 Improper drum address for store of C; 

a white switch up and down will cal- 
culate first available address and 
store C. 

FF005 location 10K Overflow on element of C. 

FF006 location 109 B sum check failure; white switch 

to try again. 

FF007 location 18F Master sum check failure; reread. 

FF008 location 0J0 Row of C read incorrectly from drum; 

white switch to try again. 
NOTE A type T parameter tape can be converted to a type D tape 

merely by raising the white switch up and down at stop 

2ko6o. 
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